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K. W. Wagner's treatment of the distribution of relaxa- 
tion times in dielectrics is reviewed ; the effect of the density 
of distribution upon the frequency variation of the dielectric 
constant and dielectric loss factor is discussed; and a 
graphical method of evaluating the constants of Wagner’s 
equation from experimental data is described. These 
constants have been evaluated for a number of typical 
dielectrics, and the frequency variations of the dielectric 
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constant and dielectric loss factor as computed from 
Wagner's equations and from the simple equations for a 
single relaxation time are compared with the observed 
behavior. The quantitative correlation obtained between 
Wagner’s equations and the experimental data investigated 
is viewed as a confirmation of Wagner's theory of a 
statistical distribution of relaxation times in dielectrics. 


INTRODUCTION 


UMEROUS investigators have found that 

the well-known simple absorption equations 
for the time variation of the reversible absorption 
current and for the frequency variation of the 
dielectric constant ¢’, and the dielectric loss 
factor ¢’, are not in accord with their experi- 
mental data. It is generally agreed that the 
failure of these simple absorption equations to 
explain the experimental data is to be ascribed to 
the polydisperse nature of the systems investi- 
gated. The form of the absorption equations is 
governed by the form of the relaxation function 
in the fundamental equations of dielectric ab- 
sorption. The simple absorption equations corre- 
spond to a relaxation function consisting of a 
single negative exponential term. Although the 


basic theories of dielectric absorption all lead to a 


relaxation function of this form, they regard a 
dielectric as being a monodisperse system, that is, 
a system in which the factors giving rise to the 
absorption are nonvariant within the medium. 
Such a system is characterized by a single 
relaxation time. Many dielectrics, however, are 
polydisperse systems, that is, the factors giving 
rise to the absorption vary within the medium 
so that the dielectric is characterized not by a 
single relaxation time but by what may be 
termed a distribution of relaxation times. 

The concept of a distribution of relaxation 
times was first introduced by von Schweidler' in 
1907. He represented the relaxation function as a 
summation of a large number of negative ex- 


* Presented at the eighth annual meeting of the Com- . 
mittee on Electrical Insulation of the National Research 
Council at Pittsfield, Mass., Oct. 17, 1935. 

1E. von Schweidler, Ann. d. Physik 24, 711 (1907). 
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ponential terms, there being one term for each 
different relaxation time. On the assumption that 
the relaxation times vary continuously from 0 to 
x, he then expressed the relaxation function as a 
definite integral involving a distribution function. 
Instead of developing a distribution function 
from theoretical considerations, he worked back- 
wards and determined mathematically the distri- 
bution function that would give the empirical 
relaxation function, 


&(t)=Bt", 


first proposed by J. Hopkinson.? 

In 1913, K. W. Wagner*® proposed that the 
distribution of relaxation times in dielectrics is 
governed by a probability function and derived 
the corresponding absorption equations for the 
time variation of the reversible absorption cur- 
rent and the frequency variation of ¢’ and ¢”’. In 
1914, he obtained quantitative agreement be- 
tween the theoretical time variation of the 
reversible absorption current and that actually 
observed in various‘samples of balata.* However, 
the validity of the equations for the frequency 
variation of e’ and e’’ apparently has never been 
investigated. In 1926, Karapetoff® proposed a 
distribution function of the form F(a) = Ha™—e-"« 
in which H, m and n are constants independent 


of a, and a@ is the reciprocal of the relaxation 


time. However, the integrals resulting from the 
introduction of this distribution function into the 


equations for the frequency variation of ¢e’ and e” 


2 J. Hopkinson, Phil. Trans. 166, 489 (1876); Phil. Mag. 
314 (1876); Phil. Trans. 167, 599 (1877); Proc. Roy. 
a 25, 496 (1876). 
3K. W. Wagner, Ann. d. Physik 40, 817 (1913). 
*K. W. Wagner, Archiv. f. Elektrotechnik 3, 83 (1914). 
> V. Karapetoff, J. A.I.E.E. 45, 236 (1926). 
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lead to infinite series when m is not an integer and 
hence very little simplification, if any, results. 
The applicability of this distribution function 
has not been investigated. 

In 1932, Whitehead and Bajfos® obtained 
excellent agreement between a charge and dis- 
charge—time relationship determined experi- 
mentally on a specimen of dry paper at 60°C and 
1500 volts and values computed from a relaxation 
function consisting of three negative exponentials 
in which the six constants were evaluated by the 
method of successive residuals as outlined by J. 
Lipka.’ Although this procedure suffices for the 
predetermination of the a—c characteristics of 
dielectrics from suitable charge and discharge 
curves, as applied by Whitehead and Bajos, it 
cannot be viewed as a confirmation of von 
Schweidler’s theory but rather as a means of 
obtaining an empirical equation by the choice of 
a sufficient number of constants. Wagner’s equa- 
tions on the other hand involve but four con- 
stants all of which have physical significance. 
The only constant not appearing in the simple 
absorption equations is b which governs the 
density of distribution of the relaxation times. 
Hence, quantitative agreement between Wagner’s 
equations and experimental data may be viewed 
not only as a confirmation of von Schweidler’s 
theory concerning a distribution of relaxation 
times but also as_a confirmation of Wagner’s 
proposal that this distribution is governed by a 
probability function. 

Recently several experimenters*™ investi- 
gating the mechanism of dielectric absorption in 
various dielectrics have found that the simple 
equations for the frequency variation of e’ and ¢’’ 
are inadequate to explain their experimental 
results which indicate broader regions of ab- 
sorption and lower dielectric loss factor maxima 


6 J. B. Whitehead and A. Bajfios, Trans. A.I.E.E. 51, 
392 (1932). 

7 J. Lipka, Graphical and Mechanical Computation (John 
Wiley & Sons, 1918). 

8H. H. Race, J. Phys. Chem. 36, 1928 (1932). 

*D. W. Kitchin, J. Ind. and Eng. Chem. 24, 549 (1932). 

1 A. H. White and S. O. Morgan, (a) Physics 2, 313 
(1932); (b) J. Frank. Inst. 216, 635 (1933). 
1935) L. Oncley and J. W. Williams, Phys. Rev. 43, 341 

12H. Rieche, Zeits. f. Physik 95, 158 (1935). 

13G. M. L. Sommerman, Rev. Sci. Inst. 5, 341 (1934); 
J. Frank. Inst. 219, 433 (1935). 

14 A. L. Ferguson, L. O. Case and G. H. Evans, J. Chem. 
Phys. 3, 285 (1935). 
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than can be accounted for by the simple ab- 
sorption equations. It is generally agreed that 
these discrepancies between experiment and 
theory are to be attributed to the polydisperse 
nature of the systems investigated. In this con- 
nection, Race* has presented an illustration in 
which he considers a hypothetical liquid dielectric 
composed of equal numbers of polar molecules of 
five different sizes. The change in the measured e’ 
would represent the sum of the contributions of 
molecules of all sizes. The loss factors, however, 
would have their maxima at different frequencies 
since the relaxation time depends upon the 
particle size. Therefore, the sum of the loss curves 
would have a maximum ordinate (e,,’’) con- 
siderably less than would have been obtained if 
the molecules had been all of one size. For the 
assumed case, the maximum ordinate of the 
sum of the five loss curves is only about 62 
percent of the maximum ordinate for the curve 
for which all the particles are assumed to be of 
one size only. The resulting ratio is nearly the 
same as that which he found between the ob- 
served and calculated maximum values of ¢’’. 

Although these authors*-" have pointed out 
that their data indicate polydisperse systems, no 
mention has been made as to the probable 
distribution of relaxation times therein nor has 
the validity of Wagner’s equations for the fre- 
quency variation of ¢’ and e” for a statistical 
distribution of relaxation times been investigated. 
The object of this paper is, therefore, to present a 
method of examining experimental data on 
anomalous dispersion and dielectric loss on the 
basis:‘of Wagner's hypothesis and to point out the 
quantitative agreement obtained between the 
experimental and theoretical curves for the fre- 
quency variation of ¢e’ and e” for various typical 
dielectrics. 


THEORETICAL DEVELOPMENT 


In order to comprehend fully the significance 
of Wagner’s equations, it is desirable to review 
briefly. the development of the fundamental 
equations of dielectric absorption. The time 
variation of the reversible absorption current, 1, 
under continuous potential is characterized by 
the expression, 


14=C,,E®(t), (1) 
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where C,, denotes the geometric capacitance; E, 
the applied voltage, and #(t), the relaxation 
function. Eq. (1) embodies Curie’s laws which 
state that the reversible absorption current is 
directly proportional to the applied voltage and 
to the dimensions of the condenser. If these laws 
are obeyed, Hopkinson’s principle of super- 
position may be applied for the determination of 
the reversible absorption current under any con- 
tinuous variation of the applied voltage E(t). 
According to von Schweidler, 7,, measured at any 
time, t,,, then becomes 


'm dE(t) 


—?)dt, (2) 
dt 


ta= ln 


where ¢ denotes the time at which a given change 
in voltage occurred. t,,—¢ then gives the elapsed 
time between the given change in voltage and the 
measurement of i,. Denoting this elapsed time 


by 4, 
dE(tm—u) 
— 


(3) 


For a sinusoidal voltage, 
dE(tm—U) 
du 


= — Ew cos w(tm—U), 


where w= 22 X frequency. Eq. (3) then becomes 


f COS w(tm— u)P(u)du. 
0 


But cos = COS wlm COS wU+SiN wl» SiN wu. 
Hence, 


c0s wl m f cos wu du 
0 


+sin wlm f “@(u) sin wu (4) 


From this equation, the dielectric loss factor is 
found to be 


cof sin wu du 


0 


(5) 


and the increase in dielectric constant due to 
dielectric absorption is found to be 
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f cos wu du. (6) 
0 


Eqs. (1), (5) and (6) may be considered the 
fundamental equations of dielectric absorption 
since they are perfectly general and apply re- 
gardless of the mechanism of dielectric absorption 
or form of the relaxation function provided that 
Curie’s laws are obeyed. Thus, if the relaxation 
function, #(u), is known, the time variation of 
the reversible absorption current and the fre- 
quency variation of the dielectric constant and 
dielectric loss factor may be determined. 

In the theories of Pellat'® and of Maxwell’ 
both as originally proposed for a_ two-layer 
dielectric and as later extended by Wagner"? for 
a disperse system, the relaxation function is a 
negative exponential of the form 


= (7) 


where k is the absorption constant and 7 is the 
relaxation time. Substituting (7) in Eqs. (1), (5) 
and (6), 


iy 
T 


=— sin wi dt, 
T 


0 
cos wt ar| 
0 


From formulae 506 and 507 of Pierce’s Table of 
Integrals, 


(la) 


= 
=e,.[ 1+k/(1+?7?) ]. 


Any mechanism of dielectric absorption in 
which the relaxation function is a simple negative 
exponential leads to equations of the form of 
(1a), (5a) and (6a). Debye’s theory'® of polar 
molecules likewise leads to equations of the form 
of (5a) and (6a) although no specific relaxation 
function is given. However, it seems reasonable 


(Sa) 


and (6a) 


18H. Pellat, Ann. Chim. Phys. 18, 150 (1899). 

16 J. C. Maxwell, Electricity and Magnetism, Vol. 1 
(Oxford, 1892), p. 328. 

17K. W. Wagner, Archiv. f. Elektrotechnik 2, 371 (1914). 

18 P. Debye, Polar Molecules (Chemical Catalog Com- 
pany, 1929). 
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to conclude that inasmuch as this theory leads 
to equations of the same form as (5a) and (6a), 
the relaxation function, if it were explicitly 
defined, would be a simple negative exponential. 

Although the various mechanisms of dielectric 
absorption all lead to a simple exponential 
relaxation function, such a simple function is not 
in accord with experimental data. In order to 
account for the disagreement between experiment 
and theory, von Schweidler introduced the con- 
cept of a distribution of relaxation times by 


extending Pellat’s theory. He regarded the. 


charged particles in a dielectric as being of two 
types, the motion of one being oscillatory and the 
motion of the other being highly damped and 
hence aperiodic. Upon the sudden application of 
an external field, the displacement of the parti- 
_ Cles of the oscillatory type is instantaneous and 
constitutes the normal displacement whereas the 
displacement of the aperiodic particles is a 
function of time, 


D,(t) = Dok(A —e-"'") 


and constitutes the so-called ‘‘viscous displace- 
ment.”’ D,(t) is the viscous displacement at any 
time t, Do is the normal displacement, & is the 
ratio of the total viscous displacement to the 
normal displacement and 7 is the relaxation time 
of the aperiodic particles. Pellat’s theory implies 
only one kind of these aperiodic particles all 
with the same relaxation time. Von Schweidler’s 
extension assumes a large number of different 
kinds of particles having different relaxation 
times. The relaxation function then consists of a 
summation of exponential terms, there being one 
term for each different relaxation time so that 


nk, 
&(t) => — ent, 


1 


k, is proportional to the number of aperiodic 
particles per unit volume having a relaxation 
time 7,. If there are an infinite number of such 
terms, with relaxation times varying continu- 
ously from 0 to «, the summation may then be 
replaced by a definite integral, 


© k(r 


T 
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in which k(r) may be regarded as the distribution 
function of the relaxation times. 
Substituting Eq. (8) in Eqs. (1), (5) and (6) 


k(r) 
i= CoE 
0 
k(r) 
arf e~'/* sin wt dt, (1b) 
0 0 T 


© k(r) 
and f arf 
0 T 


(k(r))/r is independent of ¢, hence 


© k(r)wrdr 
of (5b) 


k(r)dr 
and + | (6b) 
o 1+??? 


Although von Schweidler derived equations of 
the form of (1b), (5b) and (6b) involving a 
distribution function, he proposed no definite 
function for the distribution of relaxation times. 
Instead, he merely used the concept of a distri- 
bution of relaxation times as justification for 
using the empirical relaxation function 


&(t) = Bi- (9) 


cos wt at} 


in deriving equations for the time variation of the 
absorption current and the frequency variation of 
e’ and ¢”. In fact, he equated the theoretical 
relaxation function containing the distribution 
function to the empirical relaxation function and 
determined what the distribution function must 
be in order that the theoretical and empirical 
relaxation functions be equivalent. He found that 
the distribution function must be 


k(r) =const. 


However, it is well known that the empirical 
relaxation function of Eq. (9) is but an approxi- 
mation formula which applies only in certain 
instances over a limited time interval. The 
introduction of an appropriate distribution func- 
tion in Eq. (8) should lead to a more general 
relaxation function. Thus, K. W. Wagner* in 
1913 proposed that the distribution of relaxation 
times is governed by a probability function 


6 
4 
12 
0.6 
068 
o3 
04 a5 AR 
+4 0. 
0.2 
10°6 10-4 10°2 102 104 


Fic. 1. The effect of the distribution constant b upon the 
density of distribution of the relaxation times according to 


the relation k(r) In 


analogous to the one proposed by Wiechert'® for 
the explanation of the phenomenon of elastic 
after-effect in solids. The distribution function 


proposed by Wagner is 
k(r)dr = (kb) “dz, (10) 


where z=In 7/79. 


This law states that the logarithms of the 
relaxation times of the various terms are grouped 
about the logarithm of the most prominent 
relaxation time ro, the density of distribution 
being governed by the distribution constant } as 
illustrated in Fig. 1. When b is =2, the grouping 
of relaxation times about the most prominent 
value ro is very dense. As b decreases from 2 to 0, 
the density of distribution broadens rapidly 
becoming infinitely broad as 6-0. 

Introducing this probability function into 
Eqs. (8), (1b), (Sb) and (6b) 


V 
C.,,Ekb 
(1c) 
kb wo cosh 2b?z9u 
=——e~¥ 20° | e~ 2 (5c) 
cosh u 


19 EF, Wiechert, Wied. Ann. 50, 335 (1893). 

*° A review of Wagner’s paper appears in “A Critical 
Résumé of Recent Work on Dielectrics’’ by L. Hartshorn, 
published in J.I.E.E. 64, 1152 (1926). In this review, cos 
is used in these equations instead of cosh presumably due 
to an incorrect translation of the German notation. 
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and 


kb 
0 


cosh (26?z9—1)u 
x au|» (6c) 
cosh u 
where Zo=In wo, 
z=In 1/79, 
and u=2+29=In wr. 


These are the equations derived by Wagner for 
the relaxation function, the time variation of the 
reversible absorption current, and the frequency 
variation of e’’ and ¢’, respectively. The definite 
integrals of these equations cannot be evaluated 
in closed form and hence must be evaluated 
numerically by approximation methods. Wagner”! 
has prepared a table of values for the integral of 
Eq. (1c) for values of 6 and ¢/ro ranging from 
0 to 1 and 0 to 25, respectively, and has also 
indicated a method* of evaluating the integrals 
of Eqs. (5c) and (6c) for values of b2 1. 

The development of the Wagner equation for 
the time variation of the reversible absorption 
current has been included in this paper for the 
sake of completeness. Hereafter, the discussion 
will be limited to the frequency variation of ¢”’ 
and ¢’ as given by Eqs. (5c) and (6c). It is of 


interest to note, however, that the same four 4 
constants appear in all three equations. Hence, %, 
if these constants can be evaluated from experi-¢ , 
mental data on the time variation of the re- ¥ 


versible absorption current, the corresponding 
frequency variation of e’ and ¢’ may be de- 
termined and vice versa. 


CONSEQUENCES OF WAGNER'S EQUATIONS 


The theoretical values of e¢’/(€,k) and 


(e’—e,,)/(€..k) as computed from a numerical 


evaluation of the integrals of Eqs. (5c) and (6c) 
for various values of 6 (0.08 to 2.0) and 
Zo (—6 to +6) are presented in Tables I and II 
and their frequency dependence (in terms of Zo) is 
shown graphically in Figs. 2 and 3, respectively. 
When b= ~, Wagner's equations transform into 
Eqs. (5a) and (6a) for a single relaxation time. In 


1K. W. Wagner, E.T.Z. 45, 1279 (1913). 
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this case, the frequency band in which dielectric 
absorption occurs is approximately two decades 
and e,»’’ is equal to 0.5¢,,k. However, as pointed 
out in the introduction, the frequency band in 
which dielectric absorption is observed in typical 
dielectrics is frequently much greater than two 
decades and e,,”’ is frequently much less than 
0.5¢,,k. That both of these apparent discrepancies 
can be accounted for by Wagner’s equations can 
readily be seen by an examination of Figs. 2 and 3 
which show clearly the effect of a distribution of 
relaxation times upon the frequency variation of 
e”’ and e’, respectively. One effect of a distribution 
of relaxation times is to increase the frequency 
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band in which dielectric loss and anomalous 
dispersion occur, the band becoming wider the 
broader the density of distribution, that is, the 
lower the value of b. In fact, as b becomes 
vanishingly small, e’ takes the constant value 


=e,(1+0.5k) 


essentially independent of frequency. A second 
fundamental effect of a distribution of relaxation 
times is to progressively decrease the relative 
magnitude of ¢,,/’ from 0.5¢,,k for a single relaxa- 
tion time (b= ~) to zero for an infinitely broad 
distribution (b=0). Thus, when b=0.25, is 
approximately equal to 0.2e,,k compared to 0.5¢,,k 


TABLE I. Values of e’’/(€,,k) as computed from Wagner's equations. 


0.25 


0.1938 
0.1919 
0.1862 
0.1750 
0.1601 
0.1426 
0.1249 
0.08889 
0.05739 
0.03385 


54.60 
' 148.4 
0.00248 403.4 


Zo 
0 
0 
1. 
1. 
2. 
2. 
3. 
4. 
5 
6 


«= 


TABLE IT. Values of (e’—«,,)/(€,.k) as computed from Wagner's equations. 


2 


0.8 d 2.0 


— 

Qo 


0.999971 
0.99978 
0.99845 
0.986 
0.964 
0.946 
0.887 
0.788 
0.658 
0.500 
0.342 
0.212 
0.113 
0.054 
0.036 
0.014 
0.00155 
0.00022 
0.000029 


0.999992 
0.999942 
0.99957 
0.9968 
0.992 
0.976 
0.942 
0.866 
0.712 
0.500 
0.288 
0.134 
0.058 
0.024 
0.008 
0.0032 
0.00043 
0.000058 
0.000008 


b 
wro* 0.15 = 0.35 0.6 0.8 1.0 2.0 » 
1.00 1.00 0.1265 0.2525 | 0.3443 0.3882 0.420 0.475 0.5000 
0.607 1.65 0.1258 0.2468 | 0.3298 0.3660 0.389 0.426 0.4432 
0.368 2.72 0.1236 0.2312 | 0.2909 0.3101 0.318 0.328 0.3238 
0.223 4.48 0.1208 0.2073 | 0.2368 0.2388 0.234 0.227 0.2123 
= 0.135 7.39 0.1165 0.1801 0.1793 0.1706 0.160 0.142 0.1325 
0.082 12.18 0.1109 0.1509 | 0.1285 0.1128 0.102 0.086 0.0814 
D 0408 D OC 0.1052 0.1195 | 0.0886 0.0713 0.064 0.053 0.0498 
0.09115 0.0684 | 0.03521 | 0.02689 | 0.02346 | 0.0195 0.0183 
ee 0.07590 0.0343 | 0.01338 | 0.00995 | 0.00865 | 0.00710 | 0.00674 
0.06065 0.0154 | 0.00496 | 0.00366 | 0.00318 | 0.00263 | 0.00248 
b 
= wTo 0.08 0.15 0.25 0.35 0.6 | 2 
- 48) 0.752 0.894 0.978 0.9965 0.999903 0.999983 0.999994 
- 74, 0.714 0.851 0.954 0.9879 0.99932 0.99988 0.999955 
- 3 | 0.675 0.798 0.911 0.9646 0.9957 0.99909 0.99966 
_ 8 | 0.632 0.734 0.843 0.910 0.967 0.9936 0.99753 
- 0.611 0.701 0.799 0.874 0.947 0.985 0.9933 
- 0.589 0.662 0.749 0.820 0.910 0.959 0.9821 
- 3 0.566 0.623 0.692 0.750 0.839 0.906 0.9525 
_ 0.544 0.585 0.634 0.671 0.744 0.816 0.8808 
- 7 0.522 0.542 0.566 0.590 0.629 0.674 0.7391 
0.500 0.500 0.500 0.500 0.500 0.500 0.5000 
0.478 0.458 0.435 0.410 0.372 0.326 0.2689 
0.455 0.415 0.366 0.329 0.256 0.184 0.1192 
0.433 0.377 0.309 0.250 0.161 0.094 0.0475 
0.410 0.338 0.251 0.180 0.090 0.041 0.0179 
0.388 0.299 0.199 0.126 0.053 0.015 0.0067 
0.367 0.267 0.155 0.091 0.032 0.0064 0.00247 
0.325 0.202 0.089 0.0354 0.00430 0.00091 0.00034 
0.286 0.149 0.046 0.0121 0.00068 0.00012 0.000045 
.4 0.248 0.106 0.022 0.0035 0.000097 0.000017 0.000006 
| 


Fic. 2. A graphical representation of “Wagner's equation 
illustrating the effect of the distribution constant 6 upon 
the frequency variation of ¢’’/(€,,k)-(zo=In wro). 


when b= «. For very small values of like 
tends to become constant essentially independent 
of frequency. From Eq. (5c), it can be shown that 


lim ¢’’ =0.886¢,,Rd. 
b+0 


The width of the frequency band in which 
anomalous dispersion and dielectric loss are 
observed in a dielectric and the ratio of e,’’ to 
0.5¢,,k are, therefore, indicative of the density of 
distribution of the relaxation times in that 
dielectric. 

Other important consequences of Wagner's 
equations are the following : 


I. When w=0, 
+h) for all values of 6. (11) 
=(0 


II. At the frequency f,, at which e’’ is a maxi- 
mum, 


(12) 


=€,,(1+0.5) for all values of b, (13) 
=€,,k* P(d). 


III. When w= 


€_. =€ 
efor all values of 
=0 


From Eq. (11), it is seen that 


=€o—e,, for all values of 


e,.k is, therefore, a measure of the total dispersion. 
Furthermore, the area under the e’’ vs. In f curve 
for any value of b is likewise a measure of the 
total dispersion. Sommerman"™ has shown that 
this area is equal to (7/2)e,k regardless of the 
form of the distribution function k(r). 
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Graphical method of evaluating the constants of 
Wagner’s equations 


In order to compare the frequency variation 
of e’ and e”’ as given by Wagner’s equations with 
that actually observed in typical dielectrics in a 
region of dielectric absorption, the constants of 
Wagner’s equations, namely 7, €,, and k must 
first be evaluated from the experimental data. 
Frequently 7) and/or ¢,, may be determined 
directly from the experimental data but often 
all four constants must be evaluated by a 
graphical method. 

In the method employed in this investigation, 
a new variable 6 is introduced and defined as the 
ratio of ¢,”’ at a frequency f, to e,’’ at a frequency 
f. where f, is any arbitrarily chosen reference 
frequency and f, is arbitrarily chosen to be one 
decade lower than f, when f,<f,, and one decade 
higher than f, when f,>f,. Under these con- 
ditions, 6 will never be greater than one and will 
have the same value for +2,, where 


2:=In (fz/fm). (14) 


For a given fixed frequency ratio, 6 is independent 
of ¢,, and k, but is a function of 6 which governs 
the slope of the loss-frequency curve and of z, 
which determines the position of f, on the loss- 
frequency curve with respect tof». The theoretical 
dependence of 6 upon 6 for various values of z, for 


1.0, 


0.9) 


0.8) 


° 
-6 -5 -4 -3 = “1 2 3 
Zo 


__ Fic. 3. A graphical representation of Wagner’s equation 
illustrating the effect of the distribution constant 6 upon the 
frequency variation of - (zo =In wro). 
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TABLE III 


TIMES IN DIELECTRICS 


. The variation of 5 with b and zz. 


b 


—) 


0.8 1.0 


S|» 


0.3452 
0.2830 
0.2336 
0.1708 
0.1382 
0.1242 
0.0998 


0.2940 
0.2360 
0.1992 
0.1504 
0.1238 
0.1125 
0.0935 


a frequency ratio of one decade is given in Table 
III and presented in the form of reference curves 
in Fig. 4. 

When f,, can be determined directly from the 
experimental data, ) may be determined directly 
from the reference curves by computing 6 at 
any convenient value of 2:, e.g., 2-=0. When 
the measured frequency band does not include 
fm, both 6 and f, must be determined simul- 
taneously from two values of 6, say 6; and de, 
computed at reference frequencies f,, and fz. 
The values of 6 and z2,; giving 6; and the values 
of 6 and 2.2 giving 52 are then determined from 
the reference curves. However, 


In(fz1/f22) 


Hence b may be plotted as a function of 2:1. 
This plot will consist of two curves the inter- 
section of which determines } and 2,1. fm and r9 
may then be evaluated by Eqs. (14) and (12), 
respectively. A representative plot from the 
analysis on Vinsol at 98°C is shown in Fig. 5. 
Two reference frequencies and hence two curves 
are sufficient to determine } and 2,;, but a third 
reference frequency and curve is included as a 
check. 

Having evaluated 6 and 79, the appropriate 
theoretical values of ¢’’/(e,.k) and (e’ —e,,)/(e,,k) 
at the frequencies at which measurements have 
been made may now be obtained from a plot 
of the theoretical frequency variation of ¢’’/(e,,k) 
and of (e’—«,,)/(e,.k) for the appropriate values 
of b and 7». The interpolation curves of Figs. 6 
and 7 make it possible to construct such plots 
for any value of b=2. If Wagner’s equations 
explain the observed data correctly, then a 
plot of the observed dielectric loss factor, éops.’’, 


at various frequencies as a function of ¢’’/(e,,k) 
at corresponding frequencies should be a straight 
line passing through the origin with slope equal 
to ¢,,k. Similarly, a plot of the observed dielectric 
constant, €ops.’, at various frequencies as a func- 
tion of (e’—e,,)/(e,k) at corresponding fre- 
quencies should be a straight line whose slope is 
also ¢,,k and intercept on the ¢’ axis is ¢,,. The 
accuracy with which the points define straight 
lines is a measure of the accuracy with which 
Wagner’s equations explain the observed data. 
A representative plot from the analysis on Vinsol 
at 98°C is shown in Fig. 8. If ¢,, and en’ can be 
determined by observation from the experi- 
mental curve, k can be evaluated directly by 
Eq. (13). If €, is unknown but ¢€o and e,,’ are 
known, then both e,, and k can be evaluated by 
solving Eqs. (11) and (13) simultaneously. 

Having evaluated 70, and k, and €w’ 
can be calculated at any frequency from the cor- 
responding values of ¢’’/(e,,k) and (e’ —e,,)/(€.), 
respectively. Using the same values of «¢,,, k and 
To, €s’ and e,’ may be calculated at any frequency 
from Eqs. (5a) and (6a) for a single relaxation 
time. The observed and calculated values of e’’ 
and ¢’ may then be compared. In order to illus- 
trate the method described above, a step-by-step 
evaluation of the constants for Vinsol at 98°C 
is presented in Table IV. 


THE DISTRIBUTION OF RELAXATION TIMES IN 
TYPICAL DIELECTRICS 


As pointed out previously, the consensus of 
opinion of both the early and recent investi- 
gators is that such dielectrics for which the 
simple absorption equations fail to explain the 
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0.15 0.25 0.35 0.6 2.0 
0.903 0.774 0.645 0.4300 0.224 0.1980 
0.876 0.729 0.589 0.367 0.1782 0.1596 
0.854 0.688 0.535 0.3105 0.1495 0.1363 
0.817 0.612 0.441 0.2250 0.1228 0.1134 
0.778 0.548 0.371 0.1775 0.1100 0.1049 
0.742 0.494 0.311 0.1505 0.1056 0.1018 
0.673 0.387 0.233 0.1135 0.088 0.1002 
0.628 0.319 0.193 0.105 0.1000 
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Fic. 4. Reference curves for determining values of 6 and 
the corresponding values of 2, giving any value of 6 
computed from experimental data. 


experimental data are polydisperse systems. 
Wagner’s equations are based upon the hy- 
pothesis that polydisperse systems are charac- 
terized by a statistical distribution of relaxation 
times. Therefore, if Wagner's equations explain 
quantitatively the experimental data on a given 
dielectric, the hypothesis of a statistical dis- 
tribution of relaxation times is established for 
that dielectric and the density of the distribution 
may be determined. 

Before analyzing the data for various di- 
electrics on the basis of Wagner’s hypothesis, it 
is advisable to consider qualitatively the various 
factors which can give rise to a distribution of 
relaxation times in a dielectric. Although Wag- 
ner’s equations apply for any mechanism of 
dielectric absorption, the factor or factors giving 
rise to the distribution of relaxation times can 
be entirely different. Thus, according to the 
Maxwell-Wagner theory of interfacial polariza- 
tion, dielectric absorption occurs in some fre- 
quency band when a material of properties €2, o2 
is distributed in concentration p in a medium of 
properties provided €20; If component 
2 is distributed uniformly in the form of spheres 


in component 1 and if furthermore p<1, then 
the relaxation time is given by 


= (2ert €2)/(4a(20;+ o2)). (15) 


Now, if component 2 is considered as being not 
pure material 2 but rather a mixture of material 
2 in material 1, then e2: and ¢2 will depend upon 
the composition of the mixture. If there are n 
kinds of spheres having different compositions 
and if, furthermore, there is a distribution of 
composition about a most prominent value, 
there must necessarily also be a distribution of 
relaxation times. A perfectly uniform mixture 
would have a singular relaxation time but a non- 
uniform mixture would have a distribution of 
relaxation times the density of distribution being 
governed by the degree of nonuniformity. Wag- 
ner’s paper!’ should be consulted for a more 
detailed discussion of the distribution of relaxa- 
tion times in heterogeneous mixtures. 

According to the Debye theory of polar 
molecules, the relaxation time of a polar di- 
electric is 


t= ((ect 2)/(€+ 2))”, (16) 
=5/(2kT). (17) 


7’ is the relaxation time of the polar molecules or 
aggregates, k is Boltzmann’s constant, T is the 
absolute temperature and ¢ is a constant meas- 


\ 


63=0600 


where 


VINSOL TEMP 98°C 
032 


\ 


-2 
2x) 


Fic. 5. The graphical determination of the distribution 
constant b and the most probable relaxation time ro for 
Vinsol at 98°C. 
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uring the internal friction. In considering the 
factors which can give rise to a distribution of 
relaxation times in a polar dielectric, it will be 
assumed that the distribution of 7 is equivalent 
to a distribution of 7’.22 From Eq. (17) it is 
obvious that a distribution of 7’ in a polar 
dielectric corresponds to a distribution of ¢. 
According to Stokes’ law, the internal friction 
opposing the rotation of a sphere of radius a in 
a liquid of inner viscosity 7 is 


Although this relation is valid only for a spherical 
particle, it seems reasonable to conclude that, 
whatever the shape of a particle, ¢ will always 
be a function of its size and shape and of the 
internal viscosity opposing its rotation. Conse- 
quently, the distribution of ¢ in a polar dielectric 
is governed by the variations in these three 
factors. 

Colloidal, associated and pelymerized systems 
are composed of particles of various sizes and 
perhaps various shapes. The inner viscosity 
opposing the rotation of a particle is probably a 
function of its size, shape and position with 
respect to neighboring particles as pointed out 
by Sommerman' and Oncley and Williams." 
Sommerman regards the inner viscosity acting 
upon a spherical particle as being a function of 
the ratio of the particle radius to the inter- 
molecular distances of the medium. In solutions 


05; T 
ZA t.0 
03 
o2 
23.0 
° a2 04 06 08 10 12 14 16 18 20 


Fic. 6. Interpolation curves for determining the fre- 


quency variation of e’’/(e,k) for any value of b=2. 
Zo=In 


22 9 in Eq. (16) may be a function of 7’ so that the 
function for the distribution of r’ may not have the same 
form as that for the distribution of 7. However, the form 
of the distribution function is unimportant in considering 
the factors giving rise to the distribution. 
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Fic. 7. Interpolation curves for determining the fre- 
quency variation of (¢’—<«,,)/(€,k) for any value of b=2. 
Zo=In 


the intermolecular distances supposedly depend 
upon the solvent and also vary within it. The 
inner viscosity assumes the experimental or 
macroscopic value of viscosity for a large sphere 
but assumes a lower value for an ion or small 
polar molecule. Oncley and Williams found that 
a change in the macroscopic viscosity of floricin 
solutions by a factor of several thousand ap- 
parently does not change the inner frictional 
constant ¢ by any large amount. In order to 
account for this behavior, they offer the following 
explanation: ‘“‘Liquids of high viscosity may 
contain large molecules but they will also contain 
some smaller ones. Thus, a viscous mineral oil 
consists of long chain solid hydrocarbons dis- 
solved in shorter chain liquid hydrocarbons. 
In a system consisting of polar molecules dis- 
solved in such a solvent, we may expect the 
arrangement of the polar molecules to be such 
that the frictional forces would be as small as 
possible. Thus, the molecules will tend to arrange 
themselves in such a way that the inner frictional 
constant ¢ will be a minimum. If very large 
molecules are present, they may form pockets in 
which the smaller molecules would have nearly 
free rotational movement. Thus, a small inner 
frictional constant ¢ will result from any obser- 
vations designed to measure it. Very small 
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Fic. 8. The graphical determination of ¢,, and k for Vinsol 
at 98°C. 


molecules will also have a small inner frictional 


constant. Accordingly, an increase in macro-— 


scopic viscosity might either increase or decrease 
the inner frictional constant. It may also be 
expected that a variable time of relaxation r, 
or inner frictional constant ¢, might result even 
though the polar molecules were only of one 
kind, the difference being due to their movements 
with respect to solvent molecules of various 
sizes and shapes.” 

It is also conceivable that the individual 
molecules in an aggregate are capable of being 
displaced from their preferred orientation by the 
external field and that the internal friction 
opposing the displacement of any molecule is a 
function of the position of that molecule in the 
aggregate. Thus, it might be expected that the 
internal friction opposing the displacement of a 
molecule in the center of an aggregate is much 
greater than that opposing the displacement of a 
molecule on the surface of the aggregate. 

Although the above considerations are only 
qualitative, they serve to point out the factor or 
factors which govern the distribution of relaxa- 
tion times in dielectrics. Furthermore, the density 
of the distribution of relaxation times in a 
dielectric is a measure of the complexity of the 
constitution of the dielectric. Thus, a perfectly 
uniform mixture or a pure unassociated polar 
liquid in which the effect of the intermolecular 
forces is negligible should be characterized by a 


single relaxation time whereas a nonuniform 
mixture or an associated or polymerized polar 
liquid should be characterized by a distribution 
of relaxation times, the density of distribution 
becoming broader the greater the degree of non- 
uniformity of the mixture or the greater the 
degree of association or polymerization of the 
liquid. Some conception of the density of the dis- 
tribution of relaxation times in a dielectric may 
be gained from knowledge of the complexity of 
the dielectric and conversely some conception 
of the complexity of a dielectric may be gained 
from an experimental determination of the dis- 
tribution constant, provided, of course, that 
Wagner’s equations are applicable to the system. 
In order to test the validity of Wagner’s 
hypothesis of a statistical distribution of relaxa- 
tion times in polydisperse systems, the following 
materials were chosen as being representative of 
dielectrics of various degrees of complexity : 


1. Glycerol**—distilled. 

2. Vinsol**—a mixture of highly oxidized and highly 
polymerized abietic acid and terpenes. 

3. Abietic acid’*—in the form of fine crystals and reputed 
to be free from mineral matter and the pitches which are 
usually associated with rosin. 

4. Glycol phthalate resin No. 1—data of R. H. Kienle 
and H. H. Race.*5 

5. Resoglaz+30 percent Arochlor** No. 1254—Resoglaz 
is a polystyrene thermoplastic resin. Arochlor™®® is a 
mixture of chlorinated diphenyls. 


Glycerol was chosen as being representative of 
a pure but associated polar liquid. The single 
molecules of an associated liquid are mutually 
attracted by secondary valence forces and aggre- 
gate into complexes containing two or more of 
the single molecules. In general, the average 
number of molecules per complex is small so 
that the distribution of size of the aggregates 
and the distribution of inner viscosity impeding 
their rotation are probably relatively dense. 
Under these conditions, associated liquids should 
be characterized by a relatively dense distribu- 
tion of relaxation times. The results of the 
analysis on glycerol (Table V and Fig. 9) confirm 
the supposition that the distribution of relaxation 
times in glycerol is relatively dense and, further- 


23S. O. Morgan, Trans. Electrochem. Soc. 65, 185 (1934). 

24 Unpublished data of the author. 

25 R. H. Kienle and H. H. Race, Trans. Electrochem. 
Soc. 65 (1934). 
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RELAXATION TIMES IN DIELECTRICS 


more, show that Wagner’s equations explain the 
observed data more accurately than the equa- 
tions for a monodisperse system. At —31.5°C, 
the agreement between Wagner’s equations and 
the experimental data is within the experimental 
error and the computation errors inherent in the 
many interpolations necessary in computing 
the four constants of Wagner’s equations. At 
—43.6°C, however, both Wagner’s equations and 
the equations for a monodisperse system give 
values of ¢’’ which at high frequencies are lower 
than the experimental values. This departure at 
high frequencies is to be attributed to the 
asymmetry of the experimental e’’ vs. In f curve 
due either to an asymmetrical distribution of 
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relaxation times or to the presence of a second 
region of dielectric absorption at higher fre- 
quencies overlapping the first. Assuming that an 
asymmetrical distribution of relaxation times is 
manifestation of a nonuniform distribution of 
complex size, then the form of the experimental 
e’’ vs. In f curve would require the distribution 
of the complexes of smaller than average size to 
be broader than that of the complexes of greater 
than average size, whereas the reverse condition 
appears to be more probable. The fact that the 
geometric dielectric constant, as determined in 
this analysis, is greater than the refraction value 
suggests the presence of a second region of 
absorption at higher frequencies. In this case, 


TABLE IV. Procedure for evaluating the constants of Wagner's equations from data on Vinsol at 98°C. 


1. Experimental 


2. Compute 6 from experi- 
data 


mental data 


5 = 


f(Ke) 


0.804 
0.761 
0.658 
0.534 
0.395 


f(ke) 
fa=1 
fes=10 


61 = =0.818 
53 = €100/ €10” = 0.600 


” 
€ 


0.1502 | 0.500 
0.1460 | 0.398 
0.1410 | 0.355 
0.1351 | 0.309 
0.1250 | 0.265 
0.1176 | 0.225 
0.097 0.155 
0.073 0.100 


a 
o 


=o 


3. From Fig. 4, record the values of 6 and the correspond- 
ing values of 2, Zr2 and 2,3 giving 6), 52 and 43, re- 
spectively. 


6,=0.818 be =0.702 5;=0.600 


b b 


to 
~ 
aw 


7. Plot vs. (Fig. 8B), and vs. (e’ 
(Fig. 8A). 
Slope either curve =«,,k =5.37 
Intercept on ¢’ axis (curve A) =¢,,=3.50 
Hence k=1.53 
. Compute and from and (e’—e,,)/e,k, 
respectively, and e,’’ and ¢,’ from Eqs. (5a) and (6a), 
respectively, using values of e¢,,, ro and k evaluated; 
and compare with eons, and e’obs.. 


0.087 
0.096 
0.116 
0.130 
0.152 
0.174 
0.194 
0.220 


0.128 
0.140 
0.165 
0.185 
0.208 
0.243 
0.269 
0.305 


| 


on 


” 
€ obs. 


0.804 
0.761 
0.658 
0.534 
0.395 


® 271 + In (fri/fro) =222 — 1.10. 
321 =223 + In (fr1/ fra) =223 — 2.30. 


. Plot b vs. 22: and determine intersection of curves. 

See Fig. 5. 
Intersection: b=0.183; 2::=0.37. 

. Determine f,, and 70. 

221 = In (fri /fm) =0.37; fm = 690 cycles per second ; 
t9= 2.31 X seconds. 

. From Figs. (6) and (7) tabulate ¢’/(e,k) and 
(e’—«,,)/(€,.k) for various values of zo when 6=0.183. 
Plot as a function of f and pick off values at 1, 3, 10, 
30 and 100 kc. 


According to Wagner's equations. 
4 According to equations for single relaxation time. 


9. Plot the observed and calculated values of ¢’’ and ¢’ as 
a function of frequency. See Fig. 9. 


m | f(ke) f(ke) 
0.69 1 | 0.1492 | 0.464 
fF 1.88 3 | 0.1415 | 0.357 
3.09 10 | 0.1230 | 0.250 
1 | 5.98 5.10 30 | 0.1001 | 0.170 
3 | 5.41 | 8.40 100 | 0.073 | 0.105 
10 | 4.84 | 13.85 ; 
30 | 4.42 ] 37.6 
100 | 4.07 } 102.3 
Zn" 
| 1.7 
0.7 f(kc)| | | daw af” af 
| _08 1 0.801 | 2.51 5.98 | 5.99 | 5.23 
13 3 0.760 | 1.17 | 5.41 | 5.42 | 3.77 
-18 10 0.660 | 0.369 | 4.84 | 4.84 | 3.53 
| 205 30 0.537 | 0.123 | 4.42 | 4.41 | 3.50 
| 230 | 100 0.392 | 0.037 | 4.07 | 4.06 | 3.50 
5 
6 


the equations for the frequency variation of ¢’ 
and ¢”’ would contain two Wagner integrals, one 
for each region of absorption. At frequencies in 
the band in which the two regions overlap, 
é”’ will be the sum of the individual contributions 
of each region. Two regions of absorption in 
glycerol could be accounted for by assuming that 
there are two most probable sizes of the com- 
plexes, one being much larger than the other. 
The smaller most probable size might possibly 
correspond to a single molecule. The low and 
high frequency absorptions would then be due 
to the restrained rotation of the larger complexes 
and the smaller complexes or single molecules, 
respectively. The absorption constant of the low 
frequency absorption in glycerol is much greater 
than that of the high so that the number of 
large complexes is probably much greater than 
the number of small complexes or single mole- 
cules unless the moment of the large complexes 
is much greater than that of the small. Two 
regions of absorption in glycerol could also be 
accounted for if the high frequency absorption 
is regarded as being due to the restrained rotation 
of the OH groups. However, the existence of a 
high frequency absorption in glycerol is un- 
certain. The above discussion is given primarily 
to point out the usefulness of an analysis of this 
kind in interpreting dielectric data. 

Vinsol is representative of a dielectric of 
relatively high degree of polymerization. Accord- 
ing to the supplier it is a mixture of highly 
oxidized and highly polymerized abietic acid and 
terpenes. On cooling from the liquid state, 
Vinsol forms a glass and exhibits dielectric 
absorption due presumably to the restrained 
rotation of polar aggregates of abietic acid and 
some of the terpenes. Since both the size of the 
aggregates and the inner viscosity impeding 
their rotation probably vary between wide limits 
in a material of such complex constitution, 
Vinsol should be characterized by a very broad 
distribution of relaxation times. The results of 
the analysis on Vinsol is given in Table V and 
Fig. 10. As predicted, the distribution of relaxa- 
tion times was found to be very broad (b=0.18). 
At all three temperatures, Wagner’s equations 
explain the experimental data quantitatively 
whereas the simple absorption equations fail 
completely. Thus, ¢’” calculated from the equa- 
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tions for a monodisperse system for 100 kilo- 
cycles at 98°, 10 kilocycles at 110° and 1 kilocycle 
at 129°C is, respectively, 0.094, 3.45 and 0.031 
times as large as the experimental value. It is 
obvious that the application of these equations 
to polydisperse systems of such complexity is 
utterly hopeless. 

Abietic acid is representative of a dielectric 
whose constitution is more complex than that of 
glycerol but less complex than that of Vinsol. 
Since abietic acid polymerizes to some extent, 
the distribution of the na‘ factor in it is probably 
broader than in glycerol which merely associates 
but, on the other hand, its degree of polymeriza- 
tion is probably lower than that of Vinsol. 
Therefore, the density of the distribution of 
relaxation times in abietic acid should be inter- 
mediate between that in glycerol and Vinsol. 
The results of the analysis on abietic acid are 
given in Table V and Fig. 11. As predicted, the 
distribution of relaxation times was found to 
be moderately broad (b=0.388), being inter- 
mediate between that in glycerol and Vinsol. 
The agreement between Wagner’s equations and 
the experimental data is again quantitative 
whereas the agreement with the equations for a 
monodisperse system is very poor. 

Kienle and Race are of the opinion that the 
dielectric absorption which they observed in 
glycol phthalate resin No. 1 and other alkyd 
resins is of the Maxwell-Wagner type. They 
picture the probable structure of such resins as 
being that of a gel composed of a heterogeneous 
mixture of macromolecules or micells of different 
sizes with conducting components held in the 
structure. The distribution of components in 
such a system is undoubtedly nonuniform and 
hence a distribution of relaxation times is to be 
expected. The analysis of the data on glycol 
phthalate resin No. 1, given in Table V and Fig. 
12, shows that the distribution of relaxation 
times in this material is moderately broad 
(b=0.327), comparable to that in abietic acid. 
Wagner’s equations again explain the experi- 
mental data quantitatively whereas the equa- 
tions for a monodisperse system are unsatis- 
factory. The fact that the geometric dielectric 
constant is considerably greater than the re- 
fraction value is additional evidence in support 
of the view that a Maxwell-Wagner interfacial 


b 

T 

€, 

k 
b 

fi 

T 

k 
b 

f 

7 

€ 

f 

7 

€ 

f 

{ 

€ 

J 

J 

A 


RELAXATION TIMES IN DIELECTRICS 447 
TABLE V. Comparison of the observed and theoretical frequency variation of ¢’ and ¢' for various dielectrics 
| f(ke) | €obs.” tw * | | | 
A. Glycerol. Temp. —31.5°C 
b 1.16 
Sax. 135 ke 1 0.38 0.35 54.2 54.2 54.2 
To 1.18 10~® sec 3 1.3 1.28 1.05 54.2 54.2 54.2 
hn 6.7 10 4.1 4.18 3.50 53.7 53.5 53.9 
k 7.09 30 11.1 10.8 10.1 50.5 50.2 52.0 
€nok 47.5 100 19.5 19.9 22.7 35.3 35.9 37.3 
Temp. —43.6°C 
b 1.19 
Smax 11 ke 1 5.0 5.1 4.27 57.4 56.8 57.3 
To 1.45 X10~° sec 3 13.2 12.6 12.1 53.0 52.5 54.4 
ra 10.2 10 20.8 20.7 23.6 35.2 35.5 36.2 
kh 4.66 30 15.8 15.2 15.3 18.6 18.3 15.8 
€ok 47.5 100 8.1 6.0 5.2 10.1 11.2 10.8 
B. Vinsol. Temp. 98°C 
b 0.183 
Suez 0.69 ke 1 0.804 0.801 2.52 5.98 5.99 5.23 
To 2.31 K10™ sec 3 0.761 0.760 1.18 5.41 5.42 3.76 
On 3.50 10 0.658 0.660 0.371 4.84 4.84 3.52 
k 1.53 30 0.534 0.537 0.124 4.42 4.41 3.49 
€..R 5.37 100 0.395 0.392 0.037 4.07 4.06 3.49 
Temp. 110°C 
b 0.18 
Sumax 10 ke 1 0.680 0.680 0.534 7.A1 7.14 8.64 
To 1.59 10~ sec. 3 0.750 0.759 1.48 6.60 6.62 8.25 
vs 3.30 10 0.782 0.782 2.70 6.00 6.00 6.00 
k 1.63 30 0.764 0.763 1.62 5.44 5.42 3.84 
€.R 5.39 100 0.674 0.680 0.534 4.85 4.86 3.35 
Temp. 129°C 
b 0.147 
Smax 617 ke 1 0.353 0.336 0.011 7.86 7.89 8.45 
To 2.58 X 1077 sec. 3 0.435 0.442 0.030 7.61 7.61 8.45 
~ 2.23 10 0.551 0.554 0.108 7.23 7.23 8.45 
k 2.79 30 0.642 0.642 0.301 6.81 6.80 8.42 
€.oR 6.22 100 0.722 0.722 0.981 6.28 6.24 8.30 
C. Abietic Acid. Temp. 86.8°C 
b 0.388 
fax 70 ke 1 0.047 0.046 0.012 3.674 3.673 3.690 
To 2.27 X10~* sec. 3 0.092 0.089 0.036 3.635 3.634, 3.688 
a 2.842 10 0.158 0.158 0.119 3.545 3.545 3.674 
k 0.298 30 0.213 0.213 0.307 3.401 3.399 3.558 
oR 0.848 100 0.225 0.226 0.399 3.201 3.206 3.071 
D. Glycol Phthalate Resin No. 1. Temp. 30°C 
b 0.327 
fax 300 ke 1 1.00 0.22 0.04 16.5 16.6 16.7 
T0 5.30 X 1077 sec. a 0.68 0.58 0.11 16.2 16.4 16.7 
in 5.8 10 1.11 1.11 0.36 15.9 15.9 16.7 
k 1.88 30 1.76 1.77 1.08 15.0 15.0 16.6 
€.R 10.92 100 2.38 2.36 3.27 13.2 13.2 15.6 
300 2.62 2.62 5.46 11.2 11.3 11.3 
1000 2.28 2.32 3.00 - 93 9.1 6.7 
E. Resoglaz+30 percent Arochlor No. 1254. Temp. 25°C 
b 0.135 
max. 19ke 1 0.0397 0.0397 0.0210 2.73 2.73 2.85 
To 8.3 X10~-°* sec. 3 0.0438 0.0438 0.0620 2.70 2.71 2.84 
_ 2.45 10 0.0467 0.0462 0.1656 2.67 2.67 2.77 
k 0.164 30 0.0470 0.0463 0.1815 2.63 2.64 2.57 
€.R 0.402 100 0.0443 0.0446 0.0739 2.60 2.60 2.46 


* Theoretical values according to Wagner's equations. 
** Theoretical values according to simple absorption equations. 
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polarization is the cause of the absorption in 
this material. 


The system Resoglaz+30 percent Arochlor : 


(Table V and Fig. 13) is an example of a system 
in which a polar liquid characterized by a rela- 
tively dense distribution of relaxation times is 
dissolved in a resin characterized by an extremely 
broad distribution of relaxation times. Resoglaz 
is very highly polymerized styrene in which the 
particle size and the inner viscosity impeding 
the motion of the slightly polar particles probably 
vary between very wide limits. The fact that 
both the dielectric constant and dielectric loss 
factor of Resoglaz are essentially independent of 
frequency over a wide frequency band, a direct 
consequence of Wagner’s equations as } ap- 
proaches zero, constitutes experimental evidence 
of an extremely broad distribution of relaxation 
times in this material (b is probably of the order 
of 0.01). Arochlor No. 1254, a mixture of chlori- 
nated diphenyls, exhibits anomalous dispersion 
and dielectric loss presumably due to the re- 
strained rotation of polar molecules or aggre- 
gates. At 25°C, e’ attains its maximum value 
around 300 kilocycles and the distribution of 
relaxation times appears to be relatively dense 
(b=about 1). The absorption observed in the 
system Resoglaz+30 percent Arochlor No. 1254 
is undoubtedly due to the Arochlor but in this 
system at 25°C, e’’ attains its maximum value at 
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Fic. 10. The experimental and theoretical frequency vari- 
ation of ¢’ and e” for Vinsol at three temperatures. 


19 kilocycles and the distribution of relaxation 
times is very broad (b=0.135). The Arochlor 
molecules or aggregates dissolved in the Resoglaz 
appear to be acted upon by a greater and more 
variable internal friction than the molecules in 
the bulk material. 


CONCLUSIONS 


The results of these analyses show definitely 
that the simple absorption equations derived for 
a monodisperse system are not applicable to 
dielectrics of a polydisperse nature. ‘Although 
these simple equations give fair agreement with 
the experimental data on dielectrics of relatively 
simple constitution such as associated liquids, 
the agreement becomes increasingly poor as the 
complexity of the dielectric increases. For di- 
electrics of even moderate complexity, these 
equations are worthless. Wagner’s equations, on 
the other hand, appear to be applicable to poly- 
disperse systems of all degrees of complexity. 
The quantitative agreement obtained between 
Wagner’s equations and the experimental data 
on the various dielectrics investigated with the 
possible exception of glycerol is indeed striking 
and is viewed as a confirmation of Wagner’s 
hypothesis of a statistical distribution of relaxa- 
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tion times in these polydisperse systems. If an 
analysis of a polydisperse system on the basis of 
Wagner’s hypothesis discloses that Wagner’s 
equations are not applicable to the system, the 
indications are that the distribution function 
governing the distribution of relaxation times has 
a different form than the one proposed by 
Wagner or that there are two or more overlapping 
regions of absorption or that there are other 
disturbing factors which must be taken into 
consideration. 

Important information concerning the die- 
lectric behavior of a polydisperse system may be 
obtained by analyzing, on the basis of Wagner’s 
hypothesis, experimental data on the frequency 
variation of ¢«’ and e” over a frequency band 
which need include only a suitable portion of a 
region of absorption. This analysis will show 
whether or not Wagner’s equations are applicable 
to the system. If they are, a statistical distri- 
bution of relaxation times in the system is 
indicated and the complete curves for the fre- 
quency variation of «’ and ¢’ and the corre- 
sponding curve for the time variation of the 
reversible absorption current can be calculated, 
provided of course that there are no other 
regions of absorption in the material. The evalu- 
ation of the constants of these equations gives at 
once the static and geometric dielectric constants, 
the most probable relaxation time, the distri- 
bution constant, the frequency at which e’’ is a 
maximum and the absorption constant of the 
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Fic. 12. The experimental and theoretical frequency 
variation of ¢ and e” for glycol phthalate resin No. 1 
at 30°C. 


material, all of which are important properties 
of a dielectric. 

The distribution constant is a rough index of 
the complexity of a dielectric. It may be used as 
a qualitative measure of the degree of association 
or polymerization of a dielectric of the polar 
type and of the degree of uniformity of the 
distribution of the various components of a 
heterogeneous mixture of the Maxwell-Wagner 
type. An estimate of the complexity of a dielectric 
may therefore be obtained from an experimental 
determination of its distribution constant. An 
investigation of the relationship between the 
distribution constant and the degree of associ- 
ation or polymerization of a dielectric would be 
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of interest as it might throw some light on the 
general problem of association and might provide 
a useful tool for controling the degree of 
polymerization in the manufacture of certain 
resins. 

The inherent difficulties in ascertaining whether 
an observed absorption is of the Debye type or of 
the Maxwell-Wagner type is well known. Un- 
fortunately, an analysis of experimental data on 
the basis of Wagner’s hypothesis appears to be of 
little value in differentiating between the two 
mechanisms of dielectric absorption. Although 
the factors causing a distribution of relaxation 
times in the two mechanisms are entirely differ- 
ent, a very dense or a very broad distribution is 
theoretically possible in either. In some cases, 
however, the magnitude of ¢,, may offer a clue as 
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to the nature of the absorption. Thus, in systems 
in which e¢, is considerably greater than the 
refraction value, the absorption is more apt to be 
of the Maxwell-Wagner type, unless there is 
evidence of a second region of absorption at 
higher frequencies as in glycerol, whereas in 
systems in which ¢,, is approximately equal to the 
refraction value, the absorption is more apt to be 
of the Debye type. 

The author is greatly indebted to Dr. L. A. 
MacColl of Bell Telephone Laboratories for the 
numerical evaluations of the complex integrals of 
Wagner’s equations and other mathematical 
assistance and to Dr. S. O. Morgan, also of Bell 
Telephone Laboratories, for his interest and 
helpful suggestions in connection with this in- 
vestigation. 
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(1). Introduction: Description of fixed-focus and variable-focus electrostatic image tube; scope 
of paper. (2). Focusing properties: Measurement of the variation of image distances and 


magnification with object distance and ratio of applied voltages for both types of tubes; calcula- 
tion of potential distributions and electron paths; comparison of experimental and theoretical 
results. (3) Aberrations: Classification of aberrations; measurement of tangential and sagittal 
image surfaces for fixed-focus tube; calculation of axial (chromatic and spherical) aberrations; 
calculation of field aberrations; comparison between measurements and calculation; reduction 


of field aberrations by curving the cathode. 


1. INTRODUCTION 


N a recent paper, Drs. Zworykin and Morton! 

have described several electron optical sys- 
tems for producing, on a fluorescent screen, 
slightly magnified or unmagnified images of 
extended objects, such as photosensitive surfaces 
on which light images are projected. The combi- 
nation of such a photocathode, a suitable electron 
lens system and a fluorescent screen for observing 
the electron image, we shall call an “image tube.”’ 
In the image tubes considered in the paper 
referred to above, the focusing of the electrons 
was accomplished by purely electrostatic means. 
The simplest (‘‘fixed-focus’’) type of tube (Fig. 
1, a) consisted of two narrowly separated cylin- 
ders of equal diameter at cathode and screen 
‘potential, respectively, the juncture of the two 
being described as the “‘lens.’’ By replacing the 
cathode cylinder by a system of rings at uni- 
formly increasing or decreasing potential (Fig. 
1, b), simulating a resistive cylinder between 
whose ends a voltage is applied, it was found 
possible to obtain exact focusing for a given 
object and image distance by varying the ratio 
of the voltages applied across the system of rings 
(Ve) and that between the anode and cathode 
(V,). Finally, it was found that, in either case, 
giving the cathode a curvature such that its 
center of curvature falls into the lens (Fig. 1, c), 
greatly improves the quality of the image, 
especially in the peripheral portions. 

In the present paper, a more detailed experi- 
mental and theoretical study of the focusing 
properties and image defects of the above image 


1V. K. Zworykin and G. A. Morton, J. Opt. Soc. Am. 26, 
181 (1936). 
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tubes will be presented. It is hoped that such a 
study may help to outline the scope of these 
image tubes by pointing out the natural limits 
imposed by the focusing characteristics and 
aberrations, and suggesting in what directions 
future improvement must be sought. 


2. FOCUSING PROPERTIES 


In order to measure the variation of image 
distance v with cbject distance u, both measured 
from the ‘‘lens,” in the fixed-focus type of tube, 
the tube shown in Fig. 2 was constructed. The 


a u Vv 
b u v 
Ve 
c u v 
CATHODE LENS SCREEN 


Fic. 1. Schematic diagram of image tubes. 
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Fic. 2. Experimental tube with movable cathode 
and screen. 


cathode and screen are movable and can be 
shifted by tilting the tube. Catches are provided 
which hold the cathode in place when the tube is 
tilted in one position, and the screen in place 
when it is rotated about its axis through 180°. 
This allows the two elements to be moved 
independently. The lens-to-image distance v, for 
which the center of the image appeared sharp, as 
well as the corresponding magnification m, were 
measured for a set of cathode-to-lens distances u. 

Fig. 3 shows a corresponding modification of 
the focusing ring type tube shown in Fig. 1, b. 
Here, the screen is made movable, and the 
cathode fixed at a distance of 1.8 lens radii from 
the lens, the image distance and magnification 
being measured for different values of V2/V;, the 
ratio of the voltage applied between the cathode 
and outermost ring, and that applied between 
the cathode and anode. 

To obtain similar information by calculation, 
it is necessary first to find the variation of the 
refractive index n, of the medium, which is 
conveniently given by the square root of the 
potential, if we choose the latter to be zero at 
points where the velocity of the electrons van- 
ishes (e.g., at the cathode if we neglect the initial 
velocities).? Then, with the aid of Fermat’s law, 


5 fnds =0, (1) 


ds being an element of path of the ray considered, 
and the variation applying to paths with fixed 
endpoints, the paths actually taken by electrons 
emitted from the cathode and, hence, the position 
and magnification of the image can be determined. 

The variation of potential within the tube is 
governed by Laplace’s equation and conditions 
imposed at the boundaries, i.e., the potentials 
applied to the cylinders and the cathode. We 
shall, for convenience, imagine the anode cylinder 


*See e.g., E. Briiche and O. Scherzer, Geometrische 
Elektronenoptik (Springer, 1934), p. 3. 
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Va 


Fic. 3. Variable focus tube with movable screen. 


to be extended to infinity. This introduces no 
appreciable error, as electrons which have entered 
the anode cylinder, owing to their high velocities, 
are but slightly affected by any ‘‘end-effects,”’ 
including among these the presence of the fluo- 
rescent screen. 7 

Cylindrical coordinates, with their origin at 
the center of the cathode, are most convenient. z 
will denote the distance from this point along the 
axis of symmetry, and 7 the separation of any 
point from this same axis. The radius of the 
focusing cylinders (‘‘lens radius”) will be taken 
as unit of length. The substitution 


x(r, 2) = F(z) -G(r) (2) 
in Laplace’s equation 
(3) 


reduces its solution to that of the two ordinary 
differential equations. 


(1/F)(d?F/d2*) = —k’, (4a) 
(1/rG) (d/dr)(rdG/dr) =k’, (4b) 


k* being a separation parameter, with the general 
solutions 
F=ae**+be—*:, (Sa) 


G=cJo(tkr) +d No(ikr). (Sb) 


The solution for the potential must have the 
form of a linear combination of products of such 
functions with the same value of k, which may 
be any value in the complex plane. The require- 
ment that the potential remain finite as z=, 
eliminates terms with complex k. Furthermore, 
the condition that it be finite on the axis, 
eliminates those containing the Neumann func- 
tion No. Finally, the fact that the potential 
vanishes (in the case of the flat cathode) all over 
the cathode plane, limits the trigonometric 
functions to be considered to sines. The solution 
for the potential thus becomes 
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g(r, 2) = f b,.Jo(tkr) sin kedk. (6) 


The coefficient b, can be evaluated by substi- 
tuting the values of g(r, z) for r=1 by the usual 
methods used for evaluating Fourier coefficients. 


For the fixed-focus tube, assuming the gap. 


between the cathode and anode cylinders to be 
negligible, the potential along the axis becomes 


2 cos ku sin kz 
= ¢(0, z) =— V; dk, (7) 
k Jo(tk) 


the first two derivatives being given corre- 
spondingly by 


2 cos kz 
#'(z) = cos ku———dk, (7’) 
0 J (tk) 
—2 sin kz 
©” (z) =—— “Wik cos ku dk. (7") 
v0 J 


Fig. 4 shows the axial potential distribution of | 


the fixed-focus type of tube with u=1.8, together 
with its first four derivatives. 

Similarly, for the focusing ring type tube, 
replacing the set of rings by a cylinder of 
uniformly increasing or decreasing potential, the 
result for the potential along the axis is 


= Ve sin ku 


+(Vi- 


cos ku\ sin kz 
) (8) 


k Io(ik) 


For any given object distance u, it is thus 
possible, by numerical or graphical quadrature of 
the above integrals, to find the potential along 
the axis for different values of z. Furthermore, 
the potential off the axis is given in terms of the 
axial potential and its derivatives by the relation® 


2 


Replacing, temporarily, the cylindrical co- 
ordinates r, 3 by Cartesian coordinates x, y 
(r?=x*+y*), Fermat’s law can be written 


3 Reference 2, p. 66, Eq. (16). 


f F(x, y, x’, y', 8)ds=6 f 8) 


(10) 


where a prime denotes differentiation with respect 
to z. The corresponding Euler equations for the 
electron paths, 


(d/dz) F, — F,=0, (11a) 


(d/dz) F, — F,=0, (11b) 
thus become 


2 + 2 
2¢ i+y’ 


(12a) 

+ -——, 
or r dz r2¢ 


| 


x’ y” 
y= -—y (14 ) 
it+x” 2» 1+x” 
(12b) 


or 1+x” r2eor 
When the rays considered are in meridian 
planes only, attention may be confined to the 


first of these equations, and x may be replaced 
throughout by 7, y being set equal to zero. If, 


Fic. 4. Axial potential for fixed-focus tube and its first 
four derivatives. 


G. A. MORTON AND E. G. RAMBERG 


Fic. 5. Electron paths. 


finally, the expression (9) is substituted for g we 
may drop, for “paraxial rays,’’ i.e., such that 
deviate infinitesimally from the axis of the 
system, all terms involving higher powers than 
the first of r and r’ yielding : 


= — (13) 


If the radial distance from the center of the 
cathode to the starting point of an electron ray, 
r(0), and its initial slope, r’(0), are given, this 
equation, by a process of numerical integration, 
enables us to find the distance of the ray from 
the axis, r(z), for any value of z. 

To locate the image position, we shall consider 
the path of an electron leaving the center of the 
cathode with an infinitesimal initial radial 
velocity. The determination of the path is simpli- 
fied by introducing, as dependent variable, 


b=1/22—1/r(dr/dz). (14) 


which reduces Eq. (13) to a first-order equation. 
The term —(dr/dz)/r is the convergence of the 
ray at the point z, i.e., the reciprocal of the 
segment of the axis measured from the abscissa 
of the point considered to the intersection with 
the tangent to the ray ; —1/2z2, on the other hand, 
is the initial convergence of the ray, which 
describes a parabola of the form r? = pz(p =a con- 
stant) in the uniform field immediately in front 
of the cathode. The convergence variable, 6, is 
thus chosen so as to vanish for z=0, where the 
convergence itself becomes negatively infinite. . 

Substituting (14) in (13), the equation govern- 
ing b becomes 


db 
dz z 2 
This equation, again, cannot be integrated 


analytically but must be evaluated numerically. 
The value of } is determined, by a step-by-step 


22 


2h 22 


(15) 


- process of integration, for some particular value 


of z, far enough to the right of the lens so that the 
ray path is essentially a straight line. From the 
convergence of the ray at this point the i image 


position is located. 


To determine the magnification, as illustrated 
in Fig. 5, the point of intersection of the image 
plane, as determined above, and the path of an 
electron leaving the cathode, with zero initial 
velocity, at a radial distance ‘‘O”’ from the center 
of the cathode, is located. The radial distance ‘‘I”’ 
of this intersection, divided by the distance “‘O,” 
gives the magnification. To calculate this path | 
Eq. (13) must be used. 

A comparison of the results of the measure- 
ments and calculations for the fixed and variable 
focus type tubes is given by Figs. 6 to 9. The 
relatively large deviations for very small image 
distances, v, may be ascribed to the interference 
of the screen with the field distribution within the 
lens. The closeness with which the magnification 
is given by v/2u is striking for both types of 
tubes. 


3. ABERRATIONS 


The aberrations which detract from the quality 
of the electron image observed on the fluorescent 
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Fic. 6. Focusing relation (fixed focus). 
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Fic. 7. Magnification (fixed focus). 


screen can conveniently be separated into two 
classes; axial aberrations and field aberrations. 
The former, chromatic and spherical aberration 
in the customary nomenclature, affect the quality 
of the image on the axis, as well as in other parts 
of the field, while the latter, in first line astig- 
matism and curvature of field, become noticeable 
only in the outer portions of the image. It will be 
seen that these constitute the most serious 
defects of the systems here studied. 

An experimental investigation of all the image 
defects of electron optical systems has been 
carried out for a point source by Diels and Knoll,‘ 
todemonstrate their identity with the aberrations 
encountered in light optics. H. Johannson’ has 
made a special study of the aberrations of the 
electron microscope. On the theoretical side, 
Glaser® has developed the theory of electron 
optical aberrations very completely, though not 
in a form in which it is directly applicable to our 
case, where the electrons enter the refracting 
field with infinitesimal velocities. Henneberg and 
Recknagel’? have, making suitable simplifying 
assumptions, given a theoretical treatment of the 


938) Diels and M. Knoll, Zeits. f. Tech. Physik 16, 617 
1 a 
( 934) Johannson, Ann. d. Physik 18, 385 (1933); 21, 274 
1934). 

®W. Glaser, Zeits. f. Physik 80, 451; 81, 647; 83, 104 
(1933); 97, 177 (1935); Ann. d. Physik 18, 557 (1933). 

7W. Henneberg and A. Recknagel, Zeits. f. Tech. 
Physik 16, 230 (1935). 
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Fic. 8. Focusing relation (variable focus). 


chromatic (and spherical) aberrations of various 
types of image tubes. 

It was not practical, in the present case, to 
carry out measurements of the axial aberrations, 
as these are not of sufficient magnitude to be 
readily determined except for so low accelerating 
voltages that no bright images could be obtained 
on the screen. For the measurement of the field 
aberrations, curvature and astigmatism, the 


os Ss 


Fic. 9. Magnification (variable focus). 
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Object: u=1.7R 


Image: v=3.46R 


R=1.5" 


v=2.18R v=1.80R v=1.57R 


Fic. 10b. 


Fic. 10. Sections through image space of flat cathode tube for various distances from lens 2. 


procedure was the following. In the tube of Fig. 2 
the object distance was kept fixed at the value 
u=1.8 and pictures were taken of the image on 
the fluorescent screen at various distances from 
the lens; such a series of photographs is shown in 
Fig. 10. The object consisted of a grid of fine 
vertical lines superposed on an array of black 
squares. The point along a vertical line through 
the optical axis of the system at which the fine 
lines appear sharpest, indicates the location of a 
sagittal, that along a horizontal line through the 
axis where they appear sharpest, that of a 
tangential image point, just as in ordinary optics. 
Having located the sagittal and tangential image 
points in every plane, it is possible to construct 
also the sagittal and tangential image surfaces. 
In order to calculate the axial aberrations, we 
consider pencils of rays leaving the center of the 
cathode with a fixed initial kinetic energy eA® 


under different angles, e.g., 0°, 45° and 90°, and 
find their points of convergence on the axis. 

To treat the first case, 0=0, we replace ® in 
Eq. (13) by #+A®% and then introduce as new 
dependent variable (which vanishes for z=0) : 


1 1dr 


A®b A® A® 
#'(0)} Le’(o)L 


r dz 
(16) 
If by is a solution of Eq. (15) corresponding to 
the case of zero initial velocity and Ab=b—Dy is 
small up to the point where ray paths become 
substantially straight (¢=2,), we can write for 
Ab(2z;): 


0 z’ 


(18) 


where A= —(1/2+6/26) +2bp, 
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A® dbo 


dz 3 42? ’’(0)z 


To find the solutions for rays with finite initial 
lateral velocity components, we introduce the 
variable Ar=r—ry where ro is the solution of the 
paraxial Eq. (13) for an equal initial lateral 
velocity and integrate the equation 


Ar’ Ar— Ar’ 
4(6+A®) 2(@+A®) 
Pry 
4(>+A0)L 86” 


[ 
ro Tor ( =) 
2(@+ae)L  \4(o+ae) 40’ 


24’ 
This is obtained from Eq. (12a) if terms involving 
third powers of r and 7’ are retained. The proper 
initial conditions follow in every case from the 
form of the parabolic path described by the 
electron in the uniform accelerating field im- 
mediately in front of the cathode. 
For the case of the fixed-focus tube with unity 
magnification and with Aé/V,;=0.0001 (e.g., 
Ab=0.4 volt, V:=4000 volts), the axial devia- 
tions from the focal point corresponding to 


Fic. 11. Schematic diagram of lens structure and sagittal and tangential pencils. 


infinitesimal initial velocities are found to be 
0.44, 0.24, and —0.03 lens radii for 09 =0°, 45°, 
and 90°, respectively. As the solution of Eq. (13) 
for our particular case indicates that the vertex 
angle of a beam, with initial lateral velocity 
[ 2(e/m)A®, }!, is equal to 2.17 (A®,/V;)! the di- 
ameter A of the circle of least diffusion is found to 
be about 0.003 lens radii, or 0.004’"(=0.1 mm), 
increasing approximately in proportion to the in- 
itial energy of the electrons. The formula given by 
Henneberg and Recknagel’ yields A = 2mA®/’ (0) 
=0.004 lens radii or 0.006’(=0.15 mm). It is 
based on considering the image tube field as 
consisting of a uniform initial field terminated by 
a short lens. It is seen that the axial aberrations, 
while not important in the present case, may be 
of significance when the field before the cathode 
is reduced (as when a very small image is desired), 
or in cases when the initial velocities are on the 
average substantially larger and the accelerating 
voltage is reduced. 

Of greatest importance are, however, the field 
aberrations, image curvature and astigmatism. 
To calculate these, the paths of a set of principal 
rays, corresponding to electrons starting from 
rest on the cathode at various distances from the 
axis, must first be obtained by integrating the 
equation 


Pry 
4b 4b 80” 


+1" (20) 
4b 40’ 26’ 


which is found by substituting Eq. (9) in Eq. 
(12a), retaining terms up to the third power in r 
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and ¢’. Pencils of rays starting from the same 
point on the cathode, with infinitesimal lateral 
velocities, will in general come to a focus only if 
they lie either in the meridional plane or a plane 
normal thereto, the focal points being termed the 
tangential and sagittal image points, respectively. 
These image points, corresponding to all possible 
principal rays, form the tangential and sagittal 
image surfaces which are the object of the 
calculation. The separation of the two surfaces 
gives a measure of the astigmatism, their curva- 
ture determines the image curvature. Fig. 11 
illustrates the tangential and sagittal pencils 
belonging to a particular principal ray, as well as 
the two image surfaces T and S. 

To obtain the tangential image point corre- 
sponding to a given principal ray ro(z), which is a 
solution of Eq. (20) for ro’(0)=0, we introduce 
the difference variable Ar=r—ro and, further- 
more, the convergence variable 


b=1/22—(1/Ar)Ar’, 


Piy ) 
) +n | 
2b 2” 
~-0)| 1 +r? 


4b 
+ +37,’ 


(21) 


which obeys the equation 


db 1 
dz 42? 


b 
(22) 


This can be integrated numerically in the same 


Ti 
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fashion as Eq. (15) giving the tangential image 
point on the principal ray. The calculation of the 
tangential image curvature on the axis can readily 
be reduced to quadrature if Ab=b—bo, where by 
is the solution of Eq. (15), is intretueed in 
Eq. (22). 

To find the sagittal image point, the second 
Euler equation, (12b), is utilized, considering the 
displacement y as infinitesimal; then Eq. (12b) 
becomes 


y"=y 
1+r,” 2¢ 


(1 +r") 
or 


Introducing 


as well as the expansion of the potential (9), a 
convergence equation analogous to (22) is ob- 


tained: 
db o” dy 

+ +n" | 
22 26 \ 2s 4h 40’ 


dz 46 4b 
b 
24’ 


(24) 


The numerical work and the obtention of the 
sagittal radius of curvature on the axis proceed as 
in the case of the tangential pencil. 

In Fig. 12, the measured image points are 
compared with the calculated image surfaces, 
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Fic. 12. Curvature of image field. S, calculated sagittal image surface; T, calculated tangential 
image surface; 0, measured sagittal image points; x, measured tangential points. 
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Fic. 13. Large diameter of ellipse of diffusion due to 
curvature of field as function of radial distance of object 
point. 


computations being made for object points on 
the cathode separated by 0.2, 0.4 and 0.6 lens 
radii from the axis, as well as for the radii of 
curvature on the axis. Vertical and horizontal 
scales are identical so that the figure gives a true 
picture of the appearance of the image surfaces. 
Fig. 13 gives the corresponding large diameters 
of the ellipses of diffusion in the plane of the 
- paraxial image point neglecting axial aberrations. 
The axial radii of curvature of the tangential and 
sagittal image surfaces are only 0.044 and 0.075 
lens radii, or 1.7 and 2.9 mm, respectively. 

It was found, as already reported in the earlier 
paper,’ that this very serious defect could be 
largely removed by giving the cathode a curva- 
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Flat Cathode 


Fic. 14. Comparison of images from tubes with flat . 
(left side) and curved (right side) cathodes. 


Curved Cathode 


ture such that its center of curvature fell, 
approximately, into the center of the lens. The 
comparison of the images given by the flat and 
the curved cathode tube, shown in Fig. 14, brings 
this out strikingly. Without, as yet, giving a 
quantitative treatment of the curved cathode 
tube, the main reason for the improvement may 
readily be pointed out. By curving the cathode in 
the above mentioned fashion, the principal rays 
are practically undeflected in the first part of 
their path, and the sagittal and tangential pencils 
of all the rays encounter practically identical 
variations of potential. As, due to the low 
velocities of the electrons in this region, their 
paths are most easily influenced in the early 
portion of their transit, it is clear that this must 
minimize the astigmatism, i.e., the difference 
between the convergence of the tangential and 
sagittal rays, and similarly reduce the curvature 
due to the similar conditions along different 
principal rays. 

In conclusion, the authors wish to express 
their thanks to Dr. V. K. Zworykin, the director 
of the Electronic Research Laboratory, and to 
the members of the laboratory staff whose 
cooperation made this work possible. 
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Directional Characteristics of Any Antenna Over a Plane Earth 
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Simple expressions for the field as a function of direction at distances of many wave-lengths 
from any antenna above a plane earth are obtained. Except near the surface of the earth the 
field varies as 1/r. In agreement with results of others it is found that on the surface of the 
earth for values of kr, (k=2x/d), large compared to unity and yet small compared to the 
(complex) dielectric constant, the vertical component of the field varies as 1/r while if the 
distance is large compared to the dielectric constant the field decreases as 1/(kr)*?. The various 
constants etc., depending on antenna form and ground characteristics, are shown to be con- 
nected with quantities previously used for the determination of radiation resistance so that the 
same numerical work used to find the total radiation from an antenna also determines its 


directional characteristics. 


N the present paper simple expressions valid at 


large distances are obtained for the intensity 
of radiation from any antenna over an arbitrary 
earth as a function of angle. This work is an 
extension of the theory reported in two recent 
papers.': * There exact formulae were given which 
enable the computation of the total power 
radiated from any antenna above a plane earth 
of arbitrary characteristics. Moreover it was 
shown that these formulae are eminently suitable 
for practical calculation. As an intermediate 
step in obtaining the total power radiated 
expressions for the fields were obtained but these 
expressions were integrals utterly impossible for 
numerical work; so little was said about them. 
We now find, however, that it is a simple matter 
to evaluate these expressions in terms of approxi- 
mations valid at large distances and so we are 
able to get the directive patterns. Thus both 
the radiation resistance and the directive charac- 
teristics of any antenna are now obtainable from 
the same mathematical formalism and such 
numerical integrations etc., as are needed to 
give either of these properties also suffice to 
give the other. 


icot /2 


{ 


where the integration is from 0 to 2/2 to 
a2/2+i and the various symbols are defined as 


= W. Hansen and J. G. Beckerley, 220 
2 W. W. Hansen and J. G. Beckerley, to appear in Dec. 
issue Proc. I.R.E. 
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To be more exact as to the results, we find 
three separate expressions for the fields at large 
distances. These expressions have different 
regions of validity and taken together cover 
most of the interesting cases. The first and lead- 
ing term in the asymptotic expression gives the 
field many wave-lengths from the origin for all 
angles. At the surface of the earth, however, 
this expression vanishes and we need the next 
term of the series. Thus we are led to examine 
the propagation of the waves at the earth’s 
surface. This has been done before by a number 
of people and the brief treatment here adds only 
a little that is new. Specifically, in addition to 
the expression mentioned above, we find two 
expansions good at the surface of the earth, one 
valid at a distance (in wave-lengths) from the 
antenna large compared to one and yet small 
compared to the square of the ratio of the wave- 
length in air to that in the ground; the other 
valid when the distance is large compared to 
both unity and this squared ratio. 

To understand the proposed method of attack 
we need the previous results. It was shown that 
the electric field is given by 


:O)+(fs, n(0’) +asfs, — 6") : 6’) 


sin 6’d6’, (1) 
follows. 
Let &,(poz : 0’) J. (kp sin 6’), (2a) 
and then Ao,=¥V Xk.&,=(1/k)¥ XAzn, 
Azn=(1/Rk)¥ X Aon, (2b) 
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where p, ¢, 2 are cylindrical coordinates with 
origin at the surface of the earth and with the z 
axis vertical. The “‘oscillator strengths” f are 
given in terms of the current distribution i(p¢z) 
as follows: 


fan(6’) =1/(2nk? sin? 6’) f 
(3) 
fan(6’) =1/ (mk? sin? 6’) f i-As,dr. 


Finally the @ are reflection coefficients which 
depend on the characteristics of the ground and 
on 6’. 


u cos 6’ —(e’u—sin? 6’)! 


ag= ’ 
u cos 6’ +(e’ u—sin? 6’)! 


(4) 


cos 6’ — (e’u—sin? 6’)! 


a3:= ’ 

cos 6’+(e’u—sin? 6’)! 
with =e+(ia/kc), (4a) 
where e¢, yw, o are the dielectric constant, perme- 
ability and conductivity, respectively, of the 
ground. 

By varying 6’ we get an infinite number of 
functions Ae, A;, all of which are solutions of 
the wave equation. Thus Eq. (1) says that the 


vector potential is a sum of such solutions with © 


the various solutions weighted by a factor g 


Z2n(0’) =fon(0’) +a2(0’)fe, 6’) ’ 
+as(O')fs, 


We may note for future reference that since 
a(r/2)= —1 we have g(7/2) =0. 

Now the exact expression (1) for E is not 
usable because, for large values of 7, the integrand 
oscillates rapidly and numerical integration is 
impossible. The problem, then, is to find an 
approximate expression for (1), valid for large 
distances from the antenna and not too hard to 
compute. We have found several methods for 
doing this. The one we will describe follows a 
scheme first introduced by Whittaker* and later 
used in connection with radiation problems by 
Weyl‘and others. The wave functions £,,(p¢z : 6’) 


(5) 


3E,. T. Whittaker, Math. Ann. 57, 333 (1902). See also, 
Whittaker and Watson, Modern Analysis, fourth edition 
(Cambridge, 1927), p. 388 and ff. 

4H. Weyl, Ann. d. Physik 60, 481 (1919). 
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from which the A are generated are considered 
as made up of an infinity of plane waves going 
in all directions with properly chosen amplitudes 
and phases. Thus the field is expressed as a sum 
of plane waves. If now we go a long distance 
from the origin in a certain direction it turns 
out that the component plane waves going in, or 
nearly in, that direction have nearly the same 
phase and add up to be important, while the 
waves going in appreciably different directions 
have practically random phases and add to 
naught. Thus to find the field in a given direction 
we need only the amplitude of the component 
plane wave going in that direction and a simple 
expression for this amplitude is easily obtained. 

To carry through the method sketched above 
we proceed as follows. First of all we avoid 
vector integrations by going back to the defini- 
tion (2b) of the A. Also we simplify the equations 


by considering for the moment only one term 
of (1). 


— 


icoot /2 
x. n(9’)E,( poz : 6’) sin 0’dé’. 
2c (6) 


The wave function £ may be built up by plane 
waves in the following manner : 


1 +r 
f (7) 
_, 
plane 
~< b 
A 
d a 
4 


Fic, 2. 


| 
| \K 
re) 
| 
Y= ZERO 
Fic. 1. 
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where 


X =z cos sin cos ¢’ 
+y sin 6’ sin ¢’. (7a) 
So we have 


1 +r /2 


sin (8) 


Now if we wish the field E at a large distance 
from the origin in a direction specified by polar 
coordinates 6, @ it is convenient to state the 
direction 6’, ¢’ in which one of the component 
plane waves is traveling by means of polar 
angles measured from an axis pointing in the 
6, direction. We therefore introduce such 
- angles, call them w, y, defined by 


cos w=cos 6 cos 6’+sin @ sin 6’ cos (¢— ¢’) (9) 
sin sin (@¢—¢’) =sin w sin : 

and illustrated in Fig. 1. Then on putting in 

expressions for x, y, 2 in terms of the spherical 

coordinates, r, 0, ¢, with polar axis coincident 

with the z axis and making the substitution (9), 

(8) becomes 


E=———vw xk, f (10) 
2c 


where »=cos w. As to the limits we first notice 
that, provided 240, it can be shown that in 
the wine integral (1) the part of the path from 


4/2 to “tie gives a contribution to the total 


that vanishes like r~*’? and is therefore in general 
negligible when kr>1. Neglecting this, then, the 
integration in (10) is over the half-sphere 
<r. 

We now perform the integration over yp after 
first distorting the path of integration as shown 
in Fig. 2. This distortion is allowable since 
g(@’) has no singularities inside the contour. 
The point of the distortion is that aleng the 
contour from a to 6 the integrand does not 
oscillate but decreases exponentially. The inte- 
gration from b to ¢ gives nothing because the 
integrand is vanishingly small and from c to d 
the integrand is zero because g(z/2)=0. More- 
over, if kr is sufficiently large the relevant 


W. HANSEN 


region of integration is correspondingly close to 
u=1 and g(6’)e'"*’ may be considered as constant 
with value corresponding to w=0, ie., 0=6’, 
¢=¢', and then the actual integration is an 
elementary one, with result 


2 eikr 
E Xk.— (11) 
2c kr 


Just how large kr must be is worth further con- 
sideration. The restrictions may be divided into 
three classes, those that depend on the method 
used, those that depend on the antenna form, 
and finally those that depend on the ground 
characteristics. For the first it is plain that we 
must have 1<kr and it is also fairly evident 
that there must be a number of waves between 
w=0 and @’=7/2; this leads to the restriction 
1/(x/2—0)*<kr. This last could perhaps be 
removed by a slightly changed derivation. Re- 
strictions of the second and third types come 
from essentially the same requirement, namely 
that the average of g(6’)e'"* over the region 
covered by the first wave on the w, y sphere 
shall not be greatly different from the value at 
w=0. Plainly this involves the ratio of the second 
derivative to the function itself. Considering the 
variation due to e‘"® it is readily verified that 
we must have n?<kr. Now in order to have an 
fn of important size it is essential that the 
antenna extend at least far enough from the 
axis to have kp~n: In other words, measuring 
all distances in wave-lengths, we must be at a 
distance greater than the square of the distance 
from the z axis to the farthest part of the antenna. 
Much the same situation is met in considering 
variations in g that are due to changes in f. 
Thus it is impossible to get rapid variations of f 
without a high antenna and this leads to a 
requirement that kr be large compared to the 
square of the maximum value of kz on the 
antenna. These last two requirements and the 
very first one too are best lumped together into 
the single statement that kr must be large com- 
pared to (kro+1)? where ro is the radius of the 
smallest sphere centered on the origin that will 
enclose the antenna.® Finally we must consider 

5 This restriction on kr is very easily understood by 
considering the problem without ground and in spherical 


coordinates. There the field involves half-order Bessel 
functions, the order of the highest important one being 
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the changes of g due to the rapid decrease of 
1+a near 6’=7/2. On computation it is found 
that this leads to ¢’/(e’u—1)*(a/2—0)<kr for 
the @ component and a similar expression for the 
@ component with ¢’ and yw interchanged. This 
last restriction keeps the forthcoming Eq. (12) 
from being valid at the surface of the ground 
and so keeps it from predicting zero field 
strength there. 

The operation VX is now performed and the 
deleted summation restored and we obtain 


(Ket("+ g3(8) sin bei"? 
c kr” 


E= 


sin (12) 


Thus we have found an expression for the 
field strength as a function of angle, this ex- 
pression being valid for large values of kr, i.e., 
many wave-lengths from the antenna. We may 
notice that if we compute the Poynting vector 
from (12) it is obviously in the radial direction 
and that on integrating over the half-sphere 
0<6<7/2 we get an expression for the energy 
radiated upward that is identical with Eg. (17) 
of reference 2. 

As remarked above Eq. (12) is invalid at z=0 
and as this is often the most interesting region 
we next look for an expression good there. 

The problem of the variation with kr of the 
field at z=0 has already been treated at length 
by many writers in the case where the trans- 
mitter is a dipole. In fact there exist fairly 
simple approximate formulae from which fields 
at all distances can be obtained. It would be in 
the highest degree surprising if the field due to a 
complex antenna varied in an appreciably 
different way at distances kr>(kro)?, though 
just what constant factor should be used is 
perhaps not quite plain. The analysis to be 
presented verifies this guess and also gives an 
expression for the constant to be used when the 
antenna is not an infinitesimal dipole. In one 
way, however, we have not gone as far as possible 


about / = kro. Then provided kr>>/? or kr>> (kro)? the Bessel 
functions are well approximated by single terms of Hankel 
type asymptotic expansions and these combine to give an r 
factor of the form e‘*’/r. But for smaller values of kr this is 
not possible; from a mathematical point of view because 
the asymptotic expansion is not accurate, from a physical 
point of view because the waves really are longer near the 
antenna. 
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for though a single formula valid for all kr 
probably exists we find instead two approximate 
formulas, good in different regions of kr. 

Before looking for general expressions, we may 
glance at three rather obvious special cases. 
First if «’ =y=1, i.e., no ground at all, then 
a@2=a;=0, the g are not in general zero at 
and (12) is both valid and useful, 
though the validity is not obvious without a 
slightly changed derivation. Second, let »=1 
and |e’|—>* for example, by letting the ground 
become a perfect conductor. Then, as inspection 
of (4) will show, —1, +1. This makes 
go(4/2)=0, = 2f3(2/2) so causing (12) to 
indicate that the ¢ component of E vanishes at 
the earth’s surface and that the 6 component is 
just twice what it would be without ground. 
Both of these conclusions are obviously correct, 
the first because there can be no tangential 
field at a perfect cenductor, the second by the 
elementary theory of images. Finally, the third 
limiting case of gives the same 
result as the second but with the 6 and ¢ com- 
ponents interchanged. 

We will now proceed with the problem of 
finding expressions for the field at z=0. We will 
carry out the work in detail only for A», merely 
stating the result for As. First of all the ex- 
pression (1) is, strictly speaking, valid only for z 
greater than the highest part of the antenna. 
For z=0 it should be replaced by a slightly 
different equation derived from Eq. (8a), refer- 
ence 1. This differs from (1) in that, for z=0, 
there occurs instead of f(6’) +af(m— 6’) the factor 
f(r—0')(1+a). Now it is convenient to make 
the change of variables L=sin 6’, M=cos 06’. 
On doing this the factor sin 6’d0’ becomes 
LdL/M and we have 


2c/o M 


where A: has been replaced by an approximate 
expression valid for large values of kpL. Also p 
has been replaced by r as the two are identical 
for z=0. Next we perform the addition 1+ a2 
and find 


1 
1+o.=2M. (14) 
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Thus 1+ a2 goes to zero like M at L=1 so 
canceling the 1/M in (13) and avoiding an 
infinity in the integrand. This leads us to 
define the function® 


he=fo(1+a2)/M, (15) 
with which we can rewrite (13) as 


E=k,— LAL. (16) 
2c/0 


To get a qualitative understanding of the next 
step we may profitably consider the integral as 
equal to the sum of an alternating series, terms 
of which are obtained by integrating between 
successive nodes of the integrand. Now for 
large kr the Bessel function oscillates rapidly so 
the terms are very nearly equal in absolute 
value and cancel to a high order except in those 
regions where h is changing rapidly. This argu- 
ment leads us to examine the behavior of / with 
the object of finding where it is changing rapidly 
and then to transform the integral in such a 
way as to make use of the above observation. 

Examining (15) and (14) we see that # has a 
sharp peak at L=1 and also a rise when L is 
near (e’u)'. This latter we will henceforth neglect 
because, for a number of reasons, it is un- 
important when e’y is large compared to one, and 
we will confine ourselves to this case. Considering 
the region near L=1 we see that he reaches a 
peak value of 2f2(7/2)u/(e’4u—1)) at L=1 and 
drops to half this value when 6=1-—L 
= (¢’u—1)/y?. When 4 is less than this value and 
also less than one, we find by the binomial 
expansion that /z is well approximated by 


(2y)! 
On the other hand if u is large so that 6 can be 


small compared to one and yet too large for (17) 
to hold, then 


he =2fo(3/2) | (17) 


(18) 


1 
he 


®Our he is equivalent to Sommerfeld’s 1/N’. Moreover 
for a magnetic dipole at z=0, f2=constant, in which special 
case (16) will be recognized as identical with expressions 
for essentially the same thing obtained by Sommerfeld. 
See, for example, Riemann-Weber, Vol. II, p. 565. There is 
a similar connection in the electric dipole case. 


We can now evaluate (16). If it happens that 
1<kr<p?/(e’u—1) then there will be only a 
fraction of an oscillation of the integrand in the 
region covered by (17) and a considerable 
number in that where (18) is valid. We therefore 
put (18) in (16) and, since the important region 
is near L=1, we are also allowed to approximate 
the Bessel function by its asymptotic expansion, 
provided kr>n?. Doing this and using 6 as an 
integration variable, we are led to integrals of 
the general form 


f F(6) sin kré a (19) 
n kr6é—, 


where F is a slowly varying function of 6. These 
are easily evaluated by using a theorem due to 
Bromwich,’ which says that if —1<8<1 and 
certain other conditions are fulfilled, then when 


x81 F(x) sin vxdx= v8 sin 382. (20) 


We see that our qualitative prediction that the 
value and manner of variation of h near L=1 
would be critical is verified. The actual evalua- 
tion is now an elementary matter and we easily 
find after restoring the deleted summation and 
the terms in As;, 


eikr 
E= 


(21) 


2c kr 


The restrictions on kr for the two terms are 
different and in fact mutually contradictory. 
If (n+1)*<kr<e” the first term is correct, but 
it will be impossible to also have krp? (and 
also ¢’/u>1); so the second term is not to be 
used. In practice it seems likely that the first 
term only will be of use as high values of » do 
not occur naturally. 

We may compare the result of (21) with the 
various special cases mentioned above. If e/—>~ 
the first term of (21) is valid for all values of kr 
and gives a field in the 6 direction just twice 
what would exist without ground. If »=1, 


then ¢”/(ue’—1)— but in spite of this 


7T. J. la’ Bromwich, Theory of Infinite Series, 174. 
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(21) is not expected to be valid because of the 
neglect of the peak in h at L=(e’y)}. It is easily 
verified, however, that in the limit «1 this 
neglect should make the result too high by a 
factor two as in fact it does. 

If in addition to being large compared to 
(n+1)?, kr is also large compared to y?/(e’u—1), 
then approximation (17) is the proper one to use. 
The first term of this gives nothing but the 
second term after manipulations similar to those 
described above gives 


eikr 
— 
2c (kr)? e’u—1 


the first term of which is valid for (7+1)*<kr 
and ¢”/e’u—1<kr; for the second term e’ in 
these conditions is replaced by yz. 

Taken together the expressions (21) and (22) 
describe the behavior of E at z=0. In agreement 
with previously known results they say that, if 
uw=1 and ¢’>1, the field is almost entirely 
perpendicular to the earth’s surface, that it 
varies like 1/kr for kre’ and like 1/(kr)? for 
kr>e’. Moreover, in the latter case, the field is 
about «’/kr times the strength it would have if 
the inverse distance formula (21) held so that 
at any fixed distance increasing ¢’ will increase 
the field in proportion until finally (22) is sup- 
planted by (21). In the limiting case «> the 
perpendicular field is just twice what it would be 
without ground (first term of (21)) and the 


(22) 
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tangential component is zero (second term (22)). 
Finally we may note that the restriction (7/2 —6) 
>e/((eu—1)*kr) prevents (12) from ever pre- 
dicting a lower intensity near z=0 than (21) 
and (22) give there. In fact it always gives a 
higher intensity, thus the field strength in- 
creases as one goes from the earth’s surface 
upward. 

The one point that is perhaps new and useful 
is the proof that, as one might expect, the in- 
tensity at z=0 is proportional to f(7/2), i.e., to 
the intensity that would exist if the ground were 
removed. Thus to find the intensity at z=0 
due to a given antenna it is only necessary to 
multiply the current distribution by A(p¢z : 7/2) 
and integrate. This is especially easy because 
when 6’= 7/2 all the A are independent of z and 
this greatly simplifies calculations. 

To illustrate the use of (12) we have prepared 
Figs. 3 and 4 which give, respectively, the field 
as a function of angle for a vertical dipole at 
the surface of the earth® and a vertical quarter- 
wave antenna grounded at the lower end. The 
needed f values we have taken from work done 
in connection with reference 2. We have 
material from the same source for a number of 
other directive patterns, but these do not have 
symmetry about the z axis and so are rather 
harder to represent graphically. For this reason 
they are omitted. 

8 An expression from which Fig. 3 could be obtained was 
found by Weyl, see reference 4. Further work has been 


done later, see especially M. J. O. Strutt, Ann. d. Physik 1, 
721 (1929) and 4, 1 (1930). 
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The Thermal Response of the Skin to Radiation* 
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The thermal sensitivity of skin to radiation has been 
studied with a new technique and the temperature changes 
in the skin surface have been accurately measured. The 
change in skin response of both whites and negroes to 


visible, infrared radiation, and sunlight has been measured. | 


Results show that the sensation evoked by the radiations 
depends upon the optical properties of the skin. White 
subjects were found to be more sensitive to heat than 
negroes. The effect of the exposed area upon the sensitivity 
was observed. The postulate that the magnitude of the sen- 
sation depends upon two factors, nerve impulse frequency 
and central association, was found to fit the present data. A 
formula for calculating the degree of association is proposed. 
The results show that the sensitivity to radiation increases 
as the exposed area is increased up to the point where 


the stimulus has dropped to 0.0004 gram cal./sec.-cm?. 
A radiation weaker than this will evoke no sensation re- 
gardless of the size of the area exposed. This was taken to 
be the threshold of the excitation of an end organ. The 
temperature change in the skin surface with this radiation 
rate was measured to be a rise in temperature at the rate of 
0.0008°C per second; sensation is perceived after 3 seconds 
of exposure. The relationship of the temperature change of 
the skin was studied. It was found that neither the change 
(AT) nor rate of change (AT/At) of skin temperature is 
adequate to explain the sensory effect. An hypothesis is 
proposed which attributes the sensation to a differential 
change in the temperature of the blood vessel net work 
within 0.1 mm of the skin surface and another at 0.3 mm. 


ODERN calorimetric studies have shown 
that the energy exchanges of the human 
body obey the first law of thermodynamics. Since 
the experiments of Atwater, Rosa and Benedict,! 
the maintenance of a balance between the heat 
produced within the body and that eliminated 
from the body surface has been recognized as one 
of the most important conditions for normal 
physiological function in man. The constancy of 
the level of body temperature (average 37°C) has 
long been used as an index of health although 
careful measurements? have revealed a diurnal 
variation of the internal body temperature of 
about 1°C, the temperature being about 36.5°C 
in the morning hours and increasing to a high 
level of 37.5°C in the evening. A resting tempera- 
ture of more than half a degree above the normal 
average is, however, usually indicative of a 
disturbance of a pathological nature. 

Body heat production arises as the result of the 
oxidation of food and of the exothermic chemical 
reactions in the body having to do with muscular 
work and other metabolic processes. The heat 
~ *From the Russell Sage Institute of Pathology in 
Affiliation with the New York Hospital, and the Depart- 
ment of Medicine, Cornell University Medical College, 
New York. 

' Atwater and Benedict, Carnegie Institution of Wash- 
ington, Publication No. 42 (1905); Atwater and Rosa, 


Report of Storrs Agricultural Experimental Station (1907), 
p. 212. 

T. Deighton, ‘Physical Factors in Body Temperature 
Maintenance and Heat Elimination,” Physiol. Rev. 13, 
431 (1933). 


from the various organs is distributed over the 
body by the blood stream and finally brought to 
the skin surface from which it is eliminated. The 
rate of heat production is known to vary with 
emotion,* exercise, state of health, age, size, sex, 
the amount and kind of food in the stomach, 
temperature, etc. A comparison of the metabolic 
rates is generally made under conditions defined 
as “‘basal,’’ that is with the subject awake and 


‘resting, 18 hours after the last meal. The heat 


production under these circumstances is known 
as the “basal metabolic rate’’ and amounts to 
about 38 kilogram calories per hour per square 
meter cf body surface for a normal man and 
about 6 percent less for a woman. During severe 
exercise the heat production may rise as high as 
ten times the basal rate but in normal men the 
heat production cannot be lowered below the 
basal level even though the body may be receiving 
from external sources more than sufficient energy 
to maintain body temperature and the. vital 
processes. 

It is now generally recognized that the body is 
provided with a nervous mechanism for the 
control of body temperature.‘ The location of the 
thermostatic mechanism is usually given as the 


lower portion of the third ventricle of the brain. 


3 E. F. Du Bois, Basal Metabolism in Health and Disease 
(Lea and Febiger, 1936). 

Martin, ‘‘Thermal Adjustment of Man and 
Animals to External Conditions,” Lancet 2, 561 (1930). 
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THERMAL RESPONSE OF SKIN 


This brain area responds to thermal changes in 
the skin by means of special receptors for 
temperature, and to temperature changes in the 
blood passing into the brain. It is the function of 
this mechanism to so regulate heat production 
and loss that a uniform body temperature of 
37°C is maintained. The control of heat pro- 
duction is limited to increasing the energy pro- 
duced. When the body becomes cold, spontaneous 
muscular contraction occurs in the form of 
shivering and teeth chattering, and the heat 
production may rise several times the basal value. 
The first response of man on exposure to a cold 
environment is the diversion of the blood stream 
away from the skin by constricting the peripheral 
blood vessels. This affords the body a consider- 
able layer of insulating matter (the thermal 
conductivity of flesh is the same as that of cork) 
and in effect “‘converts the integument into a suit 
of clothes.’’ Continued exposure, however, will so 
cool the blood flowing to the brain that shivering 
is brought on. If by means of an artificial 
injection of immune typhoid vaccine or the 
natural infusion of malarial parasites in disease, 
the thermostatic mechanism be “‘set’’ on a higher 
temperature plane, all of the above reactions to 
cold are brought on even at normal body temper- 
ature. Thus one has a “‘chill’’ necessarily followed 
by a fever. 

Upon exposure to a warm environment the 
skin blood vessels are dilated to facilitate the 
passage of heat to the surface. Continued ex- 
posure activates the sweat glands so that body 
heat is entirely dissipated by evaporation of 
moisture when the external temperature reaches 
36°C or 95°F. In this manner, under all circum- 
stances of internal heat production and usual 
external temperature change, the body tempera- 
ture is maintained within normal limits. The 
present investigation is concerned with the sensi- 
tivity of the part of the regulatory mechanism 
which responds to heat and thermal radiation. 

Du Bois, Hardy, and their co-workers*:have 
shown that under normal conditions, 24 percent 
of the heat lost from the body is due to the 
evaporation of moisture from the lungs and skin, 
between 55 percent and 65 percent is due to the 
infrared radiation from the surface, and the 


5 E. F. Du Bois and J. D. Hardy, Proc. Nat. Acad. Sci., 
Washington (1935). 
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remaining 15 to 20 percent to convection and 
conduction. It is evident that the major portion 
of the heat loss is dependent upon the level of 
skin temperature, a fact which would lead one to 
conclude that, even within the comfort zone 
(23°C to 27°C), a systematic change in the skin 

temperature with environmental temperature is 
brought about by means of a_ subconscious 
mechanism. Such a dependence of skin tempera- 

ture upon environmental change is well known 

and the present interest in this phenomenon is 
the manner in which this regulation occurs. 

The radiation to which man is exposed in his 

natural environment comprises a spectral range 
from 0.294 in the ultraviolet to 20 in the 
infrared. Physiologically, this range may be 
divided into four regions, namely the ultraviolet 

(A<0.36u), the visible (0.36u<A<0.8u), the 
penetrating infrared (0.84 <A <3), and the non- 

penetrating infrared (A>3,). The infrared spec- 
trum is so divided because of the intense ab-— 
sorption bands of water and organic compounds at 
about (3u) and beyond. The action of ultraviolet 
light upon the skin is largely of a photochemical 

character and this type of radiation does not 
occur naturally in sufficient amounts to affect the 
thermal mechanism of the body. The effects of 
the longer waved radiation is, insofar as is now 
known, entirely thermal and the spectral divi- 
sions are made on the basis of the difference in 
the optical properties of the skin for each region. 
The sun, which is the most important radiator in 
man’s environment, has an energy distribution at 
the earth’s surface on a clear day about as 
follows: 5 percent ultraviolet, 40 percent visible 
light, 54 percent penetrating infrared and 1 

percent nonpenetrating infrared. The usual hot 
bodies in man’s surroundings radiate most of 
their energy in the near infrared, and all objects 
including man himself radiate in the nonpene- 
trating infrared. The importance of these radi- 
ation exchanges is realized from a consideration 
of the fact that while lying under the tropical 
sun, a white man absorbs about three times as 
much heat from the sun as his basal heat 
production, and when exposed to zero weather 
loses heat at about twice his basal rate by 
radiation alone. This wide range of external 
conditions and the narrow limits of normal body 
temperature make severe demands upon a ther- 
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mostatic mechanism which is adequate for the 
control of the loss and production of body heat. 
The excellence of such a mechanism, in common 
with all thermostats, depends upon the rapidity 
of its response and the sensitivity of its detector 
of thermal changes. The present study concerns 
itself with the quantitative determination of the 
sensitivity of the human skin to thermal changes 
produced by radiation, and with a correlation of 
these data with the subjective sensation of 
warmth. 

Many experiments have been performed to 
determine the sensitivity of the skin to heat and 
cold.* The methods usually consist of placing hot 
and cold objects upon the skin and following the 
external and internal thermal changes by means 
of thermocouples buried in the tissue at different 
depths. From such experiments the nerve endings 
responsible for temperature sense have been 
identified as the Krause and Ruffini organs 
located 0.1 mm and 0.3 mm, respectively, below 
the skin surface. A satisfactory idea of the 
thermal changes responsible for sensation can 
hardly be obtained from experiments of this type 
as the proximity of the sensitive endings to the 
skin surface does not permit the measurement of 
the changes in thermal gradients by means of 
wire thermocouples. Furthermore, the presence 
of foreign material near the end organs is known 
to inactivate them. Also, the use of the thermo- 
couples for the measurement of surface tempera- 
ture is open to criticism.’ It is not surprising, 
therefore, that few definite conclusions could be 
drawn from data obtained with a thermocouple 
technique. The careful work of Bazett and his 
co-workers, using thermocouples, demonstrates 
the inadequacy of the method. Other experiments 
have been made with radiant energy to determine 
the maximum intensity which the skin can bear.® 
The function of skin pigment in relation to its 
effect on the absorbed radiation has received 
scant attention and it is obvious that such 
studies must be made before the question of skin 
sensitivity to radiation can be approached. 

The thermal effect of incident radiation de- 


*H. C. Bazett, B. McGlone, R. G. Williams and H. M. 
Lufkin, Arch. Neurol and Psychiat. 27, 489 (1932). 
7J. D. Hardy, “Radiation of Heat from the Human 
II,” J. Clin. Invest. 13, 605 (1934). 
* C. Sonne, ‘“‘The Mode of Action of the Universal Light 
Bath V,” Acta. Medica. Scand. 54, 374 (1921). 
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pends upon the amount of energy absorbed by 
the skin surface and the amount which penetrates 
into the body beyond the depth of the sensitive 
endings. The question of the reflecting power of 
the skin and its opacity to penetrating infrared 
and nonpenetrating infrared has been studied 
and it has been found that white and negro skin 
have the same optical properties in these regions.° 
That is, human skin absorbs 95 percent of all 
radiation of wave-length greater than 3yu within 
0.1 mm of the skin surface and reflects less than 
1 percent ; between 1y and 3y the reflecting power 
of white and negro skin rises to about 18 percent 
and of the absorbed energy only 1 percent reaches 
a depth of greater than 3 mm. For the visible 
spectrum the reflection and penetration of white 
and negro skin have received less study. It is 
known that the reflecting power of white skin 
for visible light is approximately 35 percent, and 
that the penetrating power of visible radiation is 
greater than that for any other radiation in the 
range under discussion.!° 


EXPERIMENTAL 


I. Dependence of sensation upon wave-length 
of stimulating radiation 

The sensitivity of the skin was investigated for 
sunlight, visible, penetrating infrared, and non- 
penetrating infrared radiation. The sources of 
radiation were obtained in the following manner. 
The sunlight was directed into the laboratory by 
means of a heliostat. The radiation was reflected 
from two first surface silvered mirrors, and its 
intensity was controlled with a diaphragm which 
varied the effective aperture. The visible light 
was obtained from a 1000-watt bulb after other 
radiations had been largely filtered out by passing 
the light through 9 cm of an 0.8 percent aqueous 
solution of copper sulphate. The spectral energy 
distributions of the artificial sources were de- 
termined with a spectrometer, and it was found 
that the visible radiation contained a band of 
near infrared between 1.0u and 1.24 amounting 
to 24 percent of the transmitted light, as shown 
in Fig. 1. Attempts to increase the relative 


® J. D. Hardy and C. Muschenheim, ‘The Radiation of 
Heat from the Human Body IV and V,”’ J. Clin. Invest. 13, 
817 (1934) and 15, 1 (1934). 

10H. Laurens, Physiological Effect of Radiant Energy, 
Chap. II. (Chemical Catalog Co., New York, 1933). 
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_ Fic. 1. Spectral energy distribution of sources of radia- 
tion used. Curve 1, visible light; curve 2, near or pene- 
trating infrared; curve 3, far or nonpenetrating infrared. 


proportion of visible radiation resulted in greatly 
decreasing the intensity of visible red rays. A 
1000-watt lamp also furnished the near infrared 
radiation after the rays had been passed through 
Corning’s “heat transmitting’”’ filter. The far 
infrared energy was obtained from a 250-watt 
electric stove, 12 cm in diameter, which was 
heated to a temperature between 300°C and 
400°C before use. 

The reflecting power of white and negro skin 
was measured for these radiations by comparing 
the diffuse reflection from the skin with that 
from a scraped block of magnesium carbonate. A 
thermopile was placed at a distance of 7 cm 
normal to the skin surface and the skin illumi- 
nated with very weak radiation from an angle 
of about 60°. The diffusely reflected light from 
the skin surface was compared with that from 
the carbonate which has been found to scatter in 
a manner similar to that of the skin. The results 
are recorded in Table I. 

These values are in good agreement with the 
work of other authors.!° 

In Fig. 2 is shown a diagram of the apparatus 
used to determine the minimum quantity of 
energy necessary to stimulate a sensation. A 
convenient sized aperture, 14.5 cm?, in the small 
screen limited a definite portion of the skin 
surface to the radiation. This aperture was so 
arranged that the test subject’s head would fit 
comfortably behind the screen when the subject 


TABLE I. Diffuse reflection from skin. 


Non- 
Visible Penetrating 
Surface Radiation | Sunlight Infrared Infrared 
MgCoO; 86% 81% 73% 4% 
White Skin 35 30 20 1 
Negro Skin | 16 17 17 1 
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LC 33.5 


Fic. 2. Schematic diagram of apparatus used to measure 
minimum detectable energy and elevation of skin tempera- 
ture during irradiation. S, source of radiation; H, subject's 
head; R2 thermopile for measuring intensity of incident 


radiation; R, thermopile for measuring skin temperature; 
LC, Leslie Cube. 


was seated. The large screen separated the 
subject from the operator and protected the 
measuring instruments from the temperature 
effects of the radiators. The source of radiation S, 
was mounted on a roll table so that the intensity 
of the stimulus could be altered by moving the 
source toward or away from the subject. A 
shutter which was operated silently was inter 
posed between the subject and the radiator. Care 
was taken at all times that the subject should not 
know,. by other means than heat perception, 
when stimulation was taking place. The ad- 
vantage of this method over those in which the 
subject himself was the judge of his own per- 
ception is obvious. 

Tests for the minimum perceptible amount of 
radiation were made in the following manner. 
The subject seated himself behind the screen 
with his head in position for some minutes in 
order to facilitate concentration. The operator 
then removed the shutter and irradiated the 
exposed skin surface with a stimulating amount 
of radiation for 3 seconds. This time of irradiation 
was chosen because it was found that if the 
stimulus could evoke a sensation it would do so 
within this period. For strong stimuli the subject 
signaled his perception within a second of the 
opening of the shutter, but for the faint sensations 
near the threshold the signal was given within 2 
to 4 seconds after the beginning of the exposure. 
By placing this restriction on the time of the 
response false signals were largely eliminated. 
The intensity of the stimulus was decreased by 
steps and the test repeated until the smallest 
intensity was determined to which the subject 
responded each time the shutter was opened. The 


J D. HARDY AND T. W. OPPEL 
Il. Minimum stimulus in g-cal./sec.-cm*. 
Nonpenetrating 
Subject Sunlight* Visible Light Penetrating Infrared Infrared 
Whites I II Ill I IL III I II III I II Ill 
Trial 1 0.0063 | 0.0075 | 0.0063 | 0.0082 | 0.0082 | 0.0088 | 0.0045 | 0.0043 | 0.0038 | 0.0024 | 0.0034 | 0.0015 
2 .0085 | .0083| .0080| .0080| .0086; .0092; .0045) .0045| .0039| .0026| .0029| .0016 
3 0080; .0085| .0091) .0081 .0048 -0030 |} .0030 
4 .0088 .0049 
Average .0081 0.0086 0.0045 0025 
Energy Absorbed .0056 .0056 .0037 .0025 
Negroes I Il I II I II I II 
Trial 1 0.0077 0.0074 0.0039 0.0048 0.0061 0.0064 0.0028 0.0019 
2 .0066 .0078 .0054 .0044 .0062 .0064 .0040 .0020 
3 .0060 .0092 .0035 .0039 .0063 0039 
4 .0084 .0042 .0064 0039 
Average .0072 .0081 0.0043 0.0063 0.0031 
Energy Absorbed .0060 .0065 0037 .0049 0031 


* Area irradiated =23.8 cm?. 


factors of concentration and fatigue were found 
to be of considerable importance so that a large 
number of measurements could not be made on 
the same subject at the same time. After these 
factors were taken into consideration two or 
three complete tests were sufficient to determine 
the minimum detectable energy. 

As soon as the operator had located the position 
of the source for the minimum sensation the 
subject removed his head and the radiometer Re, 
was placed in the aperture. R2 was provided with 
a metal flange which rested against the surface of 
the small screen in such a way that the thermopile 
had its absorbing surface in the same plane 
as that the skin had occupied. The sensitive 
element of Re was an eight junction circular 
thermopile (Coblenz type) whose absorbing 
surface was 9 mm in diameter and was blackened 
with India ink. The calibration of this instrument 
was carried out by means of a Leslie Cube. The 
voltage generated by the instrument was meas- 
ured with a Leeds and Northrup potentiometer 
(Type Ke), the sensitivity of the device being 
such that 0.000126 g-cal./sec.-cm? produced 
0.00001 volt. 


(The calibration of this thermopile as recorded is in error 
by a factor 3.2. Thus, 0.000039 gram cal./sec./cm? pro- 
duced 0.00001 volt, is the correct calibration. Although the 


energy values in the article are too great by this factor, 
interpretations of tables and curves are not changed. The 
one energy value whose absolute magnitude is of importance 
is that for the minimum detectable radiation which the 
body can perceive, which is 0.00012 gram cal. per second 
per cm? instead of 0.0004.—Authors.) 


Tests were made on twelve subjects who | 


showed moderate individual differences, but the 
ratios of the values for the different radiation 
bands were the same for all. The sensitivity 
values for three white subjects and two colored 
subjects shown.in Table II are representative. It 
will be seen from the table that for the white 
subjects the important difference is a gradual 
increase in sensitivity with increasing wave- 
length. This increase in effectiveness of the 
radiation parallels the opacity of the skin to the 
energy. When the white skin is painted with 
India ink the sensitivity to all wave-length bands 
is the same and equal to that of nonpenetrating 
infrared. 

The negro subjects, although less sensitive to 
thermal stimuli, showed the same changes in the 
infrared regions as the whites, and it is only in 
the visible range that any significant effect due 
to pigment is seen. This difference in the visible 
spectrum is not sufficient, however, to make the 
negro appreciably more sensitive to sun light 
than the white man. ; 
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II. Dependence of sensation upon the surface 
area irradiated 

The actual values of the stimuli given in Table 
II are not of themselves of any great interest 
since the area exposed to the radiation was chosen 
for convenience only. It has long been known 
that the sensation produced by a heat stimulus of 
fixed intensity increases with the area stimulated, 
but the matter has not received any systematic 
study due to a lack of a satisfactory technique. 
Radiation furnishes an ideal mode of stimulation 
for this investigation for the following reasons: 
(1) the strength of the stimulus can be readily 
changed and accurately measured; (2) the size 
of the area stimulated is limited only by the size 
of the experimental subject and can be easily 
measured; (3) the time of application of the 
stimulus can be measured and controlled ; (4) the 
measurements of the stimuli are objective; (5) 
the skin temperature changes associated with any 
stimulus strength can be determined (see Section 


III). 


TABLE III. Minimum stimuli in g-cal./sec.-cm? for various 
sized areas on the body surface. 


Area Minimum Stimuli Skin Temperature Rise 
Subject | Subject Aver- Rate of 3 See- 
Location Size I I] age Rise onds 
Forehead 0.95em?) 0.0268 0.0294 
.0278 .0204 0.0268 |0.060°C /sece.| 0.18°C 
3.46 .0088 0091 
.0082 d 
|0.021 0.063 
7.08 0057 .0057 
.0060 
0054 0.012 0.036 
10.0 
.0038 .0039 .0039 0.009 0.027 
14.5 0026 .0034 
.0029 - 
.0030 0031 |0.007 0.021 
23.8 .0020 
0018 |0.004 0.012 
40.0 0013 .0014 
0015 .0014 |0.003 0.009 
Entire Face | 197 .00051 .00076 
.00076 .00063 .00066 |0.001+- 0.004 
Face and I 1680 .0006 .0007 
Chest II 1840 .0005 .0008 
.0006 .0006 \0.001+- 0.004 
Back to I 1680} .0005 
Waist II 1840 0005 .0009 .00065 |0.001+- 0.004 
Anterior I 5440 .0004 0006 
Body 11 5590 | .0005 
Surface .0004 .0005 |0.001— 0.003 — 
Posterior I 5440| .0004 .0004 
Body II 5590 .0005 .00045 |0.001— 0.003 — 
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Fic. 3. Change in the energy required to evoke sensation 
with the area of the forehead exposed. 


Of the three variables to be measured, namely, 
the magnitude of the sensation, stimulus strength, 
and surface area exposed, the first was eliminated 
as far as possible by keeping the sensation con- 
stant, i.e., at the minimum value. Thereby, the 
relationship between the stimulus strength and 
the area exposed, for a constant sensation, can be 
investigated. The experiments were done with 
nonpenetrating infrared radiation using the two 
white subjects of Section I. Far infrared was used 
in these tests because of its invisibility and the 
similarity of its stimulating qualities to other 
heat stimuli such as warm water or warm objects 
placed on the skin. 

Because of its known uniform innervation the 
forehead was first studied. Starting with an area 
of 0.95 cm? the area exposed to the radiation was 
increased by steps using circular apertures of 
larger diameters up to 40 cm’. The intensity of 
the radiation necessary to evoke a sensation for 
each size of the exposed area was measured in the 
manner described in Section I. In addition to the 
forehead, the whole face, the face and thorax, the 
anterior body surface, the back, and the entire 
posterior body surface were tested. The areas of 
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Fic. 4. Change in energy required to produce sensation 
with the size of body area exposed. 


these larger surfaces were measured by means of 
a silhouette. The results of the tests are shown in 
Table III. Although the sensation results from 
the thermal changes in the skin, these changes are 
known (Section III) to be directly proportional 
to the intensity of the radiation. Therefore, in 
this discussion the radiation rate is used as a 
quantitative estimate of the stimulus strength. 

The decrease in the amount of radiation re- 
quired to produce a sensation as the exposed area 
is increased, confirms in a quantitative way the 
qualitative observations of others. The close 
agreement between the average values of the two 
subjects probably cannot be expected for every 
pair so tested, but the tests of the minimum 
stimuli made on a series of 12 subjects showed 
that a variation of more than 15 percent could 
usually be traced to lack of concentration or 
other disturbance. The values for the results on 
the forehead and face are shown in Fig. 3. 
The points fall on a smooth curve of the type of 
an equilateral hyperbola. The smoothness of 
the curve and the fact that the location of the 
test area on the forehead does not affect the 
stimulus value would lead one to conclude that 
the end organs’ distribution over the forehead 
surface is effectively uniform. These facts also 
suggest that although the nerve ending itself is 
punctiform, it responds to stimuli over a con- 
siderable area. This is probably due to the circu- 
lation of blood in the skin area near the end 
organ. The close analogy between this arrange- 
ment and the usual physical thermopile’ with 
radiation receivers is obvious. It will be noticed 
that when an area as large as the face has been 
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exposed, a sensitivity to heat has been reached 
which is not appreciably increased even by the 
exposure of the whole body. It is seen that 
relatively smali portions of the body surface 
such as the face, arm, back, etc., are each as 
sensitive as the entire body surface. The extreme 
sensitivity of the skin to heat is evident from the 
fact that a rate of rise of skin temperature of 
less than a thousandth of a degree centigrade 
per second is perceived subjectively in three 
seconds. The whole of the data in Table III are 
represented graphically in Fig. 4. 

The shape of the curve of Fig. 3 is of the usual 
exponential type and an analysis is easier in 
terms of the logarithms of the variables. There- 
fore the logarithms of the stimuli and of the areas 
are plotted in Fig. 4. The portion of the curve 
due to the forehead is a straight line. When 
the area is increased to include the face the line 
breaks sharply and becomes almost parallel to 
the area axis. Under no circumstances was a 
stimulus of less than 0.0004 g-cal./sec.-cm? ob- 
served regardless of the size of the area exposed. 

Other portions of the body surface were 
examined in order to determine whether such a 
constant relationship as that found for the fore- 
head could be found elsewhere. Tests were made 
on the arm and back which showed that unless 
counts of the end organs were made little could 
be learned. 

From the results given in Table III it is evident 
that the sensation evoked by a stimulus is due 
to the added effects of the separate end organs, 
in much the same way that the voltage produced 
by a thermopile is increased with the number of 
junctions. In order to test further the combining 
power, or spatial summation, of the endings in 
various portions of the body the following experi- 
ments were performed. Two square apertures 
(25 cm? in area, 3 cm apart) were arranged so 
that the skin areas placed behind them could be 
stimulated singly and together. The procedure 
was to test each aperture alone and then together. 
Tests on the forehead showed that the right and 
left apertures gave identical values and together 
gave the value to be expected from Fig. 4. 
Tests were made of the backs of both hands, and 
of the forehead and the back of the right hand. 
The results are shown in Table IV. 

The sensitivity of both hands together is seen 
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to be 50 percent greater than either hand alone, 
indicating that summation is taking place. The 
sensitivity of the combination of hand and 
forehead is no greater than the forehead alone, 
showing no summation. 


III. The thermal changes in the skin accom- 
panying heat sensation 


From the foregoing it is evident that the tem- 
perature changes in the skin which accompany 
the minimal sensations are extremely small so 
that the usual procedure of following these 
changes with thermocouples either on the skin 
surface or buried in the tissue would not be 
likely to furnish reliable information. The net- 
work of blood vessels in the highly heat im- 
pervious tissue would invalidate the interpreta- 
tion of thermal gradients from thermocouple 
measurements made at points deep in com- 
parison with those at which the end organs are 
located. The presence of the thermocouple wires 
in the tissue also alters the normal heat flow to 
a great extent, especially at the point at which 
the thermocouple measures the temperature. 
For these reasons it was felt that it was not 
technically feasible to make direct measurements 
of the internal temperature changes in the layer 
of vascular tissue 0.3 mm thick which provoke 
temperature sensation. 

The current theories as to the thermal changes 
responsible for heat sensation are four in number 
and are, briefly: The Hering" hypothesis which 
assumes sensation is dependent on the absolute 


TABLE IV. Summation stimuli for various body areas. 


Backs of Both Hands 


Subject I Subject II 
Sum- Sum- 
Left | Right | Both | mation | Left | Right | Both | mation 
0.0029 | 0.0039 | 0.0020 0.0039 | 0.0041 | 0.0027 48% 
.0026 | .0025| .0018 47% 
.0018 
.0020 


Back of Right Hand and Forehead 


Fore- Sum- Fore- 


Sum- 
head |Hand}| Both | mation head 


Hand} Both | mation 


0.0016 |0.0035) 0.0016 0% 


0.0019 [0.0039 0.0019} 0% 


1 Hering, Sitzungsb. der Koeniges Akad. d. Wissensch., 


Wien 75, 101 (1877). 


Fic. 5. Radiometer device used for measurement of skin 
temperature. Vertical section of radiometer (lower left); 
Leslie Cube (right), and detail of thermopile construction 
(upper right). 


change in the skin temperature (AT) ; the Weber 
hypothesis’? supposes the rate of temperature 
change (dT /dt) to be the effector; the Ebbecke 
hypothesis" which assumes sensation to be due to 
a change in the internal thermal gradients of the 
tissues within 0.5 mm of the skin surface; the 
fourth hypothesis is a modification of Ebbecke’s 
theory proposed by Bazett and McGlone who 
postulate change in blood temperature rather 
than in the tissue. 

The present experiments were designed first 
with the idea of determining the magnitude of 
the skin temperature changes accompanying 
minimal sensations of heat, and second, by means 
of changing the internal thermal gradients with 
radiations of different penetrability, of throwing 
some light upon the questions of the thermal 
changes actually responsible for the sensation of 
warmth. 

The apparatus for measuring the skin tem- 
perature is shown in Fig. 5. The method consists 
of comparing the radiation of the skin with that 
from an experimental black body.’ By this 
means, the skin surface is not touched, the 
normal heat loss of the skin is not altered, nor 
is the skin compressed or indented with conse- 


!@E. H. Weber and R. Wagner, Handwérterbuch der 
Physiologie (Braunschweig, F. Vieweg und Sohn 1846), 
Vol. 3, p. 549. 

18 Ebbecke, Arch. f. d. ges. Physiologie 169, 395 (1917). 

4H. C. Bazett and B. McGlone, “Studies in Sensation 
III,” Arch. Neurol and Psychiat. 28, 91 (1932). 

15]. D. Hardy, “Radiation of Heat from Human Body 
I,” J. Clin. Invest. 13, 593 (1934). 
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quent vascular changes as is the case when con- 
tact thermocouples and such devices are used. 
The surface of the white human skin emits 
radiation in the neighborhood of 9u like a black 
body radiator,'* so that no emissivity correction 
is necessary. The experimental black body is a 
Leslie Cube into one side of which has been built 
a conical cavity. The cone is painted black and 
serves adequately as a black body, the energy 
values determined therefrom agreeing within 1 
percent with those obtained from a Radiation 
Standard of the National Bureau of Standards. 
During the experiments this black body was 
located near (Fig. 2) the skin surface to be 
measured and was kept at a constant tempera- 
ture (33.5°C) by means of a thermostat. The 
instrument for detecting the radiation consisted 
of an eight junction bismuth alloy thermopile, 
the sensitive elements of which were arranged in 
a circular disk 1 cm in diameter. The disk was 
blackened with bismuth black and was supported 
by a thin film of nitrocellulose. In front of the 
thermopile was arranged a silver cone 10 cm 
long in order to focus the radiation from the skin 
surface. The sensitivity of the device was such 
that a difference in temperature between the 
Leslie Cube and the skin surface of 0.01°C would 
cause a deflection of 1 mm on the galvanometer 
scale. The instrument required 7 seconds to 
reach equilibrium after being removed from one 
radiant source to another. The radiometer was 
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Fic. 6. Cooling curves of skin after irradiation showing 
extrapolation back to instant of shutting off radiation. 
Curve 1=irradiation for 15 seconds; curve 2, for 30 
seconds; curve 3, for 45 seconds; curve 4, for 60 seconds. 
a cal./sec.-cm? of nonpenetrating infra- 
r 


©]. D. Hardy, “Radiation of Heat from Human Body 
III,” J. Clin. Invest. 13, 615 (1934). 
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mounted on the vertical shaft of a pendulum so 
that it could be swung very rapidly from its 
normal position in front of the Leslie Cube to a 
position facing the skin. For obvious reasons the 
temperature changes of the skin could not be 
followed during the irradiation of the skin so 
the temperature was measured immediately be- 
fore and after irradiation. 

It was soon discovered that the temperature 
changes accompanying minimal sensations of the 
large areas were too small to be measured di- 
rectly, especially when it was considered that 
the change was always perceptible after 3 seconds 
of irradiation. For these reasons the temperature 
changes with larger intensities and long times of 
irradiation were measured and extrapolation 
resorted to for the smaller intensities and for 
the short times. 

The arrangement for controlling and measuring 
the intensity of the source of radiation was the 
same as that employed for the measurement of 
the minimum stimuli. The forehead of the subject 
was used as a test area as its temperature is 
more constant than that of any other part of 
the body surface. The radiometer R, is the 
measuring device for skin temperature and the 
Leslie Cube is shown in its relation to the skin 
surface to be measured. The head of the subject 
was protected from the heat of the cube by a 
layer of insulation. The temperature of the cube 
was read by a sensitive thermometer readable 
directly to 0.005°C. 

The subject was required to sit for about 
half an hour in the room before measurements 
were begun in order that his skin might come 
into equilibrium with the prevailing atmosphere. 
When it was found possible to obtain constant 
forehead temperature (within a few hundredths 
of a degree) the measurements were begun. 
The subject’s forehead was then irradiated at a 
known rate for a fixed time. Immediately follow- 
ing the irradiation the radiometer was swung 
over to the skin surface and the cooling curve 
to the skin followed for 120 seconds. The intervals 
chosen for irradiation were 15, 30, 45 and 60 
seconds. The type of curve obtained is shown in 
Fig. 6. 

By extrapolation the cooling curves backward 
8 seconds, the elevation of skin temperature due 
to the irradiation can be read. Such cooling 
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Fic. 7. Heating of the skin with the various radiations for constant times of irradiation. Dots =15 seconds exposure; 
triangles, 30 seconds; circles, 45 seconds; crosses, 60 seconds. No. 1 =nonpenetrating infrared; No. 2= penetrating 


infrared; No. 3=visible radiation. 


curves were plotted for various intensities of 
visible, near infrared and far infrared radiation. 
The temperature changes thus obtained were 
used to plot the heating curves of the skin for 
constant irradiation as shown in Fig. 7. 

The direct proportionality between the energy 
absorbed by the skin and the rise in temperature 
is evidence of the fact that possible vasomotor 
and environmental factors did not enter in as 
serious errors of measurement. It is felt therefore 
that extrapolation from the experimental range 
of energies to smaller energies is valid. 

The tests for sensitivity of the whole body 
were made only with nonpenetrating infrared 
radiation. Interpolation of the curves of Fig. 8 
gives the temperature values shown in Table III. 

The skin temperature elevations which occur 
when the skin absorbs the same number of 
calories of nonpenetrating infrared, penetrating 
infrared, and visible radiation is shown in Fig. 8. 
It is evident that the far infrared is the most 
efficient in raising the skin temperature, corre- 
sponding to the facts that it is most efficient in 
stimulating sensation and is totally absorbed 
within 0.1 mm of the skin surface. The thermal 
effects of the other radiations also correspond to 
their stimulating and penetrating qualities. In 
order to determine the skin temperature eleva- 
tions by the radiations which produce the same 
sensation careful measurements were made of 
minimum stimuli, using an area of 7.08 cm?. 
This area was chosen because the energy values 
are well within the experimental range of the 


temperature measurements. The results are 
shown in Table V. The average values are 
believed to be sufficiently consistent to be used 
for separating the effects due to the difference in 
penetration. The heating curves plotted from 
these energies are recorded in Fig. 9. It is evident 
that an equal sensation does not mean an equal 
rate of rise of temperature. Conversely, if the 
intensities of the three sources be adjusted so 
as to produce the same skin temperature rise, 
it will be found that the greatest sensation is 
produced by the far infrared, next by the near 
infrared, and next by the visible light. Thus, 
equal sensation does not mean an equal rise in 
skin temperature. 


DISCUSSION 


The dependence of heat sensation upon wave- 
length of the stimulating radiation is shown to 
be a matter of the optical properties of the skin 
alone, that is, of the reflectivity and the extinc- 
tion coefficient of the skin for visible, penetrating 
infrared, nonpenetrating infrared radiation.” 
All differences in the stimulating effects are 
removed by covering the skin with black paint. 
Natural pigment, however, is seen to change the 
properties of the skin for visible radiation only, 
so that no idea of pigment function can be had 
by studying the effects of such artificial pigments 
as lamp black. The effect of natural skin pigment 

17 P. C. Foster, ‘Effect of Infrared on Tissue Tempera- 


ture Gradient as Influenced by Pigment,’’ Proc. 
Exper. Biol. and Med. 33, 62 (1935). 
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Fic. 8. Skin temperature change with absorption of equal 
energy of visible light, penetrating and nonpenetrating 
infrared. 


upon the response of man to his thermal environ- 
ment seems to be negligible so that negroes are 
not appreciably more sensitive to sunlight than 
whites and they respond in exactly the same way 
as do whites to infrared radiation. Pigment 
would seem therefore to have nothing to do 
with the regulation of heat exchanges between 
man and his environment, for not only does it 
not affect the sensitivity of the body to radiant 
heat but it does not aid in any way the elimina- 
tion of heat from the body surface. Thus pigment 
affords no protection against the thermal effects 
of the natural sources of radiation and its pro- 
tective action against the specific effect of ultra- 
violet radiation is doubtful, since the layers of 
the skin superficial to the pigment absorb 
practically all of the erythemal and antirachitic 
radiations.'® Although pigment formation is un- 
doubtedly stimulated by the ultraviolet radia- 
tion, the effect of the pigment seems to lie entirely 
within the visible range. 

The straight line relationship for the logarithms 
of the values of stimulus and area for the forehead 
is so definite that it would seem to justify an 
attempt at relating these observations to sen- 
sation. That is, has a constant sensation definite 
meaning in terms of stimulus strength and 

18 W. W. Coblenz and R. Stair, ‘“‘Transmissive Prop- 
erties of Eye Protective Glasses and Other Substances,” 


Technologic Papers of the National Bureau of Standards 
No. 369 (1928). 
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TABLE V. Minimum stimuli of absorbed energy in 
g-cal./cm?/sec. 
Exposed Area = 7.08 cm?. 


Visible Near Infrared) Far Infrared 
Subject 
I II I II I II 
Trial 1 0.0101 |0.0109 |0.0074 |0.0086 |0.0055 |0.0060 
2 0112} .0114| .0074| .0052| .0052 
3 0115} .0111 .0088| .0081 | .0064} .0052 
4 | .0112| .0087 .0078| .0045} .0052 
5 .0117| .0116|) .0087| .0078| .0054)| .0052 
Average -O111 |} .0112} .0081 .0081! .0050| .0053 


number of sensitive endings exposed? Although 
many of the physiological processes taking place 
between the stimulation of a sensory ending and 
the perception of sensation have been outlined by 
others, a complete picture of the sensory mecha- 
nism must necessarily include a study of the 
subjective impressions related to these processes. 
In line with modern thought, that sensation arises 
as the result of volleys of nerve impulses con- 
ducted to the brain by afferent nerve fibers from 
stimulated end organs, the following correlation 
is made between the present stimuli data of the 
forehead and the magnitude of sensation. 

It is known that the frequency of the nervous 
discharge from an end organ increases with the 
strength of stimulation,'® provided the stimulus 
be greater than the threshold of excitation of the 
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Fic. 9. Change in skin temperature by heating with 
radiation of visible light (circles), penetrating infrared 
(dots) and nonpenetrating infrared (crosses) which will 
produce the same sensation. 


19E. D. Adrian and Y. Zotterman, “‘The Impulses 
Produced by Sensory Nerve Endings III.” J. Physiol. 61, 
465 (1926). 
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ending, a phenomenon which is rather similar to 
the reaction of a photoelectric cell when exposed 
to light. The simplest assumption in regard to the 
relationship of sensation, stimulus strength, and 
frequency of nervous discharge is that the magni- 
tude of the sensation is proportional to the 
impulse frequency reaching the brain, and the 
addition or association of the impulses into 
sensation. In the present experiment there were 
two conditions related to the sensation: (1) the 
sensation was constant in magnitude, and (2) the 
magnitude was that of the minimal sensation. 


Thus, Minimum sensation = constant 
= function (impulse frequency, association) (1) 


is the relationship to be expected, on the basis of 
the above assumptions, from the present experi- 
ments. Of the factors in Eq. (1), the impulse 
frequency will depend upon: (1) the magnitude 
of the stimulus, (2) number of end organs 
stimulated, (3) the sensitivity of the individual 
endings. Factor (1) has been studied by Hart- 
Mathews" and others. Mathew’s work 
on a cat’s muscle shows that in the range of small 
stimuli the impulse frequency is roughly pro- 
portional to the logarithm of the stimulus 
strength. In accordance with Section II, the as- 
sumption is made that the end organ distribution 
over the forehead surface is uniform. That the 
sensitivity of all the endings is approximately the 
same cannot be definitely established by the 
present experiments but is evidenced by the fact 
that the minimum stimuli for such areas as the 
face, the arm and hand, the chest, the back, etc. 
are all the same. Thus, the total number of im- 
pulses arising from a stimulated area on the fore- 
head would be proportional to log gXA, ¢ being 
the intensity of the radiation and A the area 
stimulated. That this does not represent a con- 
stant sensation is seen from Table VI, where the 
tabulated values of g¢XA are seen to increase 
with area size. The factor of central association 
or spatial summation must then be considered. 
By “spatial summation” is meant the addition in 
the brain of the separate effects of each end organ 
in the production of sensation. 

The data of the present experiment can be ac- 


20H. K. Hartline and C. H. Graham, J. Cellular and 
Comparative Physiol. 5, 277 (1932). 
21 B. H. C. Mathews, J. Physiol. 71, 64 (1931). 


curately represented by the formula (see values 
of A%-78x¢), 


log g+a log A=C (2) 


where a and ¢ are constants whose values are 0.78 
and —1.502, respectively. The constant ‘‘c’’ can 
be associated directly with the magnitude of the 
sensory response, for as pointed out above, im- 
pulse frequency is known to be proportional to 
log g. The constant ‘‘a’’ we associate with sum- 
mation of the separate end organs effects. Were 
the value of ‘a’ unity, summation would be 
perfect and a constant stimulus magnitude would 
represent a constant sensory response. The value 
of ‘“‘a’”’ has two interpretations neither of which 
can be definitely established by these experi- 
ments. Either the number of end organs per unit 
area decteases systematically” with decrease in 
stimulus intensity or summation fails progres- 
sively with smaller intensities. The evidence in 
this experiment indicates that summation is the 
more important factor, and the minimum values 
of the summation are shown in Table VI. 

Comparing Eq. (2) with the known laws of 
vision**: *4 which are, 


sensory response 
=a log ¢+6 (constant area) (3) 


sensory response 
=c log A+d (constant intensity) (4) 


shows that Eq. (3) (Weber-Fechner law) and 4 
(Granit-Harper law) apply equally well to heat 
sensation. This identity affords gratifying evi- 
dence of the fundamental similarity of the pe- 
ripheral and optical neurophysiological function. 


TABLE VI. 

Area Stimulus Percent 
(cm?) | (cal./sec.-cm?) |AX@| | Summation 

0.95 0.0268 0.025 0.026 

3.46 .0091 .032 .024 75 

7.08 .0054 .038 .025 66 
10.0 .0042 .042 .025 59 
14.5 .0031 045 025 55 
23.8 .0017 045 021 61 
40.0 .0014 056 025 44 


2 Selig Hecht, J. Gen. Physiol. 11, 255 (1927). 
(1930) Granit and P. Harper, Am. J. Physiol. 95, 211 
%G. T. Fechner, Abhandl. K. Sachs. Ges. Wissensch., 
math-phys. KI. 4, 31 (1858). 
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It would seem, therefore, that an interpretation 
of sensation can be made in terms of the stimulus 
and number of end organs. The relationship 
appears to be that the magnitude of the sensa- 
tion depends upon the total frequency of im- 
pulses reaching the brain as modified by cen- 
tral association, and an equal sensation means 
an equal number of associated impulses. No 
value of association as high as 100 percent has 
been observed, and the magnitude of this vari- 
able depends probably upon habit. The function 
of summation in these experiments is analogous 
to the electrical resistance in a thermopile circuit. 
Perfect summation would mean a physiological 
thermopile of no resistance; the falling off of 
summation with increased numbers of end organs 
is therefore the physiological simile to the 
increased thermopile resistance with the number 
of junctions. 

The difference between the threshold of sen- 
sation and the threshold of excitation for an end 
organ is clearly demonstrated by these experi- 
ments. The threshold of excitation for a single 
end organ is 0.0004 cal./sec.-cm?, whereas that 
for sensation is more than 100 times as great. 
With relatively small areas a minimum stimulus 
in the neighborhood of 0.0004 cal./sec.-cm?, is 
reached and when the area is increased beyond 
such a point no change in sensitivity is realized, 
that is, the minimum stimulus is constant. The 
sensation, however, increases as the exposed area 
is enlarged so that although the sensation is no 
longer minimal any decrease in the stimulus 
causes the sensation to disappear. Thus, the limit 
of sensitivity of the body is set not by degree of 
association of impulses but by the threshold -of 
excitation of the end organ. 

The temperature changes in the skin which 
produce sensation are seen from Section III to 
depend not only upon the heat absorbed at the 
skin surface but also upon the amount of heat 
absorbed in the various layers of skin below the 
surface. Thus, neither the actual change in 
surface temperature nor its rate of change can be 
assumed to be thermal changes responsible for 
heat sensation. Therefore, the Hering and Weber 
hypotheses may be ruled out. Body tissue itself 
‘is a very good insulator and this fact, together 
with the results given in Section III for small 
areas, which indicate that the end organs act as 


D. HARDY AND T. W. OPPEL 


though supplied with a highly conducting re- 
ceiving area, would rule out the possibility that 
the endings respond to changes in tissue tempera- 
ture alone or that such changes are usually the 
important ones. Ebbecke’s unmodified hypothesis 


can thus be said to be inadequate, although it is 


evident that his postulate in regard to a change 
in the internal thermal gradients being respon- 
sible for sensation is upheld by the present 
experiments. The hypothesis which comes nearest 
to explaining the present data is that the thermal 
receptors are really differential end organs, one 
sensitive spot located in a vascular bed superficial 
to the other, both elements being in close 
proximity to and surrounded by blood vessels. A 
differential change in the temperature of these 
blood vessel net works is assumed to be the 
necessary stimulus. A somewhat similar theory 
was proposed by Bazett and McGlone although 
their results seem to support a chemical rather 
than a thermal gradient as the stimulus. 

The above interpretation is open to the criti- 
cism that no actual measurements were made of 
the temperature changes below the skin surface. 
As we have stated, we do not believe that 
significant measurements of these changes can 
be made directly. There is ample evidence, how- 
ever, that the thermal gradients are changed 
differently by the types of radiation used in the 
above experiments and differences in the thermal 
sensitivity of the skin are found to correspond to 
these changes in thermal gradients. There is one 
phenomenon associated with temperature sen- 
sation which is particularly contradictory to any 
gradient hypothesis. Bazett found that when a 
thin fold of skin, containing end organs on only 
one side, was stimulated from the opposite side, 
a warm stimulus produced a warm sensation and 
a cold stimulus a cold sensation. On the basis of 
the gradient hypothesis paradoxical sensations 
would be expected such as are experienced when 
the skin is stimulated with a cool stimulus after 
having been subjected to a cold environment for 
some time. Since we have not repeated Bazett’s 
experiment, our results must be considered as 
only confirmatory evidence for the gradient 
theory. 

The magnitude of the temperature changes 
which produce minimal subjective sensations of 
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warmth are shown in Table III. A total rise in 
skin temperature of less than 0.003°C produced 
at a rate of 0.001°C per second is the sensitivity 
of the body as a whole and of many of its larger 
parts. Such sensitivity is evidence of the re- 
markable detector mechanism for the thermal 
regulator of the body temperature, and compares 
favorably with the most delicate of artificial 
thermoregulators. It is interesting to note the 
effect of habit upon the sensitivity of parts of the 
body. The parts of the skin surface usually most 
exposed to thermal changes are the face and 
hands and these parts although relatively small 
in area have a sensitivity almost equal to that of 
the whole body. 

The radiation sensitivity of the skin (510-4 
g-cal./sec..cm?) is not great when compared 
to that of a modern vacuum thermopile which 
can detect 10-" g-cal./sec.-cm®. It does 
compare favorably with the usual Coblenz type 
thermopile which is used in air. The thermopile 
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used in the above experiments to measure the 
radiant energy, produced about 2.5 drum divi- 
sions at the minimum stimulus. A change in the 
temperature of the walls of a room of 3°C (that 
is a 30 percent change in the normal radiation 
rate), air temperature, velocity, humidity, etc., 
being constant, will evoke a sensation in -3 
seconds. The energy exchange in this process is 
one-thirty millionth of the normal hourly radi- 
ation loss. That a sensation need not be necessary 
to this regulatory process is evident as end organ 
impulses can go to the brain without evoking a 
sensation. These impulses may have connection 
with the vasomotor pathways so that adjustment 
of internal thermal gradients of the skin can be 
made unconsciously, thereby decreasing end 
organ activity. In this way the small changes in 
environmental conditions can be automatically 
compensated for, and warning given of the larger 
changes which can be guarded against by the 
addition of clothing, etc. 
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